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X ■ Abstract: 



We present an improved calculation of -B — ^ light pseudoscalar formfactors from light-cone 
sum rules, including one-loop radiative corrections to twist-2 and twist-3 contributions, 
and leading order twist-4 corrections. The total theoretical uncertainty of our results at 
zero momentum transfer is 10 to 13% and can be improved, at least in part, by reducing 
the uncertainty of hadronic input parameters, in particular those describing the twist-2 
distribution amplitudes of the vr, K and 77. We present our results in a way which details 
the dependence of the formfactors on these parameters and facilitates the incorporation 
of future updates of their values from e.g. lattice calculations. 
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1 Introduction 

This paper aims to give a new and more precise determination of the decay formfactors 
of B mesons into hght pseudoscalar mesons, i.e. vr, K and 77. We do not include the r]' 
which is too heavy to be treated in this framework. The calculation uses the method of 
QCD sum rules on the light-cone, which in the past has been rather successfully applied 
to various problems in heavy-meson physics, cf. Refs. [1, 2, 3, 4, 5]^; an outline of the 
method will be given below. Our calculation improves on our previous papers [3, 4] by 

• including radiative corrections to twist-3 contributions to one-loop accuracy, for all 
formfactors; 

• a precisely defined method for fixing the sum rule specific parameters; 

• using updated values for input parameters; 

• a careful analysis of the uncertainties of the formfactors at zero momentum transfer; 

• a new parametrization of the dependence of the formfactors on momentum transfer, 
which is consistent with the constraints from analyticity and heavy-quark expansion; 

• detailing the dependence of formfactors on nonperturbative hadronic parameters 
describing the n, K, 77 mesons, the so-called Gegenbauer moments, which facilitates 
the incorporation of future updates of their numerical values and also allows a con- 
sistent treatment of their effect on nonleptonic decays treated in QCD factorisation. 

The motivation for this study is twofold and related to the overall aim of B physics 
to provide precision determinations of quark flavor mixing parameters in the Standard 
Model. Quark flavor mixing is governed by the unitary CKM matrix which depends on 
four parameters: three angles and one phase. The constraints from unitarity can be 
visualized by the so-called unitarity triangles (UT); the one that is relevant for B physics 
is under intense experimental study. The over-determination of the sides and angles of 
this triangle from a multitude of processes will answer the question whether there is new 
physics in flavor-changing processes and where it manifests itself. One of the sides of the 
UT is given by the ratio of CKM matrix elements iKtft/Kbl- |Kfe| is known to about 2% 
accuracy from both inclusive and exclusive b —>■ ciu transitions [7], whereas the present 
error on \Vub\ is much larger and around 15%. Its reduction requires an improvement of 
experimental statistics, which is under way at the B factories BaBar and Belle, but also 
and in particular an improvement of the theoretical prediction for associated semileptonic 
spectra and decay rates. This is the first motivation for our study of the B ^ it decay 
formfactor f^^'^, which, in conjunction with alternative calculations, in particular from 
lattice [8], will help to reduce the uncertainty from exclusive semileptonic determinations 
of \Vub\- Secondly, formfactors of general B flight meson transitions are also needed 
as ingredients in the analysis of nonleptonic two-body B decays, e.g. B -^ Kn, in the 
framework of QCD factorization [9], again with the objective to extract CKM parameters. 
One issue calling for particular attention in this context is the effect of SU(3) breaking, 
which enters both the formfactors and the K and 77 meson distribution amplitudes figuring 



^See also Ref. [6] for reviews. 



in the factorization analysis. We would like to stress here that the implementation of 
SU(3) breaking in the light-cone sum rules approach to formfactors is precisely the same 
as in QCD factorization and is encoded in the difference between vr, K and 77 distribution 
amplitudes, so that the use of formfactors calculated from light-cone sum rules together 
with the corresponding meson distribution amplitudes in factorization formulas allows a 
unified and controlled approach to the assessment of SU(3) breaking effects in nonleptonic 
B decays. 

As we shall detail below, QCD sum rules on the light-cone allow the calculation of 
formfactors in a kinematic regime where the final state meson has large energy in the rest- 
system of the decaying B, E ^ Aqcd- This is in contrast to lattice calculations which 
presently are available only for i? — »• vr and g^ > 15GeV^, due to the restriction to tt 
energies smaller than the inverse lattice spacing.^ First unquenched results are underway 
[10, 11], which, once published, will allow one to exploit the complementarity of lattice 
simulations and light-cone sum rules in more detail. 

The physics underlying B decays into light mesons at large momentum transfer can be 
understood qualitatively in the framework of hard exclusive QCD processes, pioneered by 
Brodsky and Lepage et al. [12]. The hard scale in B decays is mb and one can show that 
to leading order in l/rrif, the decay is described by two different parton configurations: 
one where all quarks have large momenta and the momentum transfer happens via the 
exchange of a hard gluon, the so-called hard-gluon exchange, and a second one where 
one quark is soft and does interact with the other partons only via soft-gluon exchange, 
the so-called soft or Feynman-mechanism. The consistent treatment of both effects in a 
framework based on factorization, i.e. the clean separation of perturbatively calculable 
hard contributions from nonperturbative "wave functions", is highly nontrivial and has 
spurred the development of SCET, an effective field theory which aims to separate the 
two relevant large mass scales m^ and ^mbAqcD in a systematic way [13]. In this ap- 
proach formfactors can indeed be split into a calculable factorizable part which roughly 
corresponds to the hard-gluon exchange contributions, and a nonfactorizable one, which 
includes the soft contributions and cannot be calculated within the SCET framework [14]. 
Predictions obtained in this approach then typically aim to eliminate the soft part and 
take the form of relations between two or more formfactors whose difference is expressed 
in terms of factorizable contributions. 

The above discussion highlights the need for a calculational method that allows nu- 
merical predictions while treating both hard and soft contributions on the same footing. 
It is precisely QCD sum rules on the light-cone (LCSRs) that accomplish this task. LC- 
SRs can be viewed as an extension of the original method of QCD sum rules devised by 
Shifman, Vainshtein and Zakharov (SVZ) [15], which was designed to determine proper- 
ties of ground-state hadrons at zero or low momentum transfer, to the regime of large 
momentum transfer. QCD sum rules combine the concepts of operator product expan- 
sion, dispersive representations of correlation functions and quark-hadron duality in an 
ingenuous way that allows the calculation of the properties of non-excited hadron-states 
with a very reasonable theoretical uncertainty. In the context of weak-decay formfac- 



^This situation may change in the future with the successful implementation of "moving NRQCD" 
[10], where the B decays while moving "backwards", which gives access to smaller values of q^ without 
increasing the discretisation error. 



tors, the basic quantity is the correlation function of the weak current and a current 
with the quantum numbers of the B meson, evaluated between the vacuum and a light 
meson. For large (negative) virtualities of these currents, the correlation function is, in 
coordinate-space, dominated by distances close to the light-cone and can be discussed 
in the framework of light-cone expansion. In contrast to the short-distance expansion 
employed by conventional QCD sum rules a la SVZ where nonperturbative effects are en- 
coded in vacuum expectation values of local operators with vacuum quantum numbers, the 
condensates, LCSRs rely on the factorization of the underlying correlation function into 
genuinely nonperturbative and universal hadron distribution amplitudes (DAs) </> which 
are convoluted with process-dependent amplitudes T. The latter are the analogues of the 
Wilson-coefficients in the short-distance expansion and can be calculated in perturbation 
theory. The light-cone expansion then reads, schematically: 



correlation function ~ ^ T^") (g) (f)^"-\ (1) 



The sum runs over contributions with increasing twist, labelled by n, which are suppressed 
by increasing powers of, roughly speaking, the virtualities of the involved currents. The 
same correlation function can, on the other hand, be written as a dispersion-relation, in 
the virtuality of the current coupling to the B meson. Equating dispersion-representation 
and the light-cone expansion, and separating the B meson contribution from that of higher 
one- and multi-particle states using quark-hadron duality, one obtains a relation for the 
formfactor describing the decay B — »■ light meson. 

Our paper is organized as follows: in Sec. 2 we define all relevant quantities, in par- 
ticular correlation functions and meson distribution amplitudes. In Sec. 3 we outline our 
calculations and derive the light-cone sum rules. In Sec. 4 we present our numerical re- 
sults and Sec. 5 contains a summary and conclusions. Detailed expressions for distribution 
amplitudes and explicit formulas for the light-cone sum rules are given in the appendices. 

2 Definitions 

The formfactors f^, f^ and fr^ which are relevant for the B ^ P transition, where P 
stands for vr, K or rj, are defined as follows:^ 

{Pip)K\B{ps)) ={(p + j>B).-^^^^^g4/W) + {^^^^ (2) 

{P{p)\J^nB{pB)) = ^ {q'{p + PB),- {ml -ml)q,}f^{q\^^), (3) 

'^ TTlB + mp 



where V^'^ = w^fjb is the standard weak current, Vj is given by V^ = 57^6 and J^ 



da^^q'^b, J^'"^ = sa^^q'^b are penguin currents. The momentum transfer is given by 



q = Pb — P and the physical range in g^ is < g^ < {rriB — mp)'^. The formfactors f^ 
and /(f are independent of the renormalization scale /i since V^ is a physical current, in 
contrast to the penguin current J°^ . Note that /+ (0) = /(f (0) which is a consequence of 

•^The following notations are frequently used in the literature: /+ = Fi and /o = Fq. 



the parainetrization chosen in Eq. (2). We assume SU(2) isospin symmetry throughout 
this work, i.e. we do not distinguish B^ -^ 7r+ and B~ -^ n^ formfactors etc. 

In the semileptonic decay B -^ irlui the formfactor f^ enters proportional to the lepton 
mass mf and hence is irrelevant for light leptons (/ = e,/i), where only /^ matters. The 
semileptonic decay can be used to determine the size of the CKM matrix element \Vub\ 
from the spectrum 

riv n^ IT/ 1 2 

where A(x) = (x + m^ — m^)^ — Arri^rn^. The formfactor f^ will be relevant in and can 
be measured from the decay B -^ ixrvr- It i^ relevant for the rare decay B -^ ttI~^1~, 
where the penguin current features in the effective Hamiltonian of the process. 

Our starting point for calculating the formfactors /^ g is the correlation function 

^,{q,PB) = zJd'xe'^-^{7r{p)\TV,{x)jl{0)\0) (5) 

= ^+(g^p|)(p + pi?)^ + ^_(g^p|)g^ , 

where Jb = ruhdi'-f^b is the interpolating field for the B meson. For the calculation of f^, 
Vfj, has to be replaced by J^. For virtualities 

ml-pl> 0(AQCD"^fe), ml-q^ > 0(Aqcd"^6), (6) 

the correlation function (5) is dominated by light-like distances and therefore accessible 
to an expansion around the light-cone. The above conditions can be understood by de- 
manding that the exponential factor in (5) vary only slowly. The light-cone expansion is 
performed by integrating out the transverse and "minus" degrees of freedom and leaving 
only the longitudinal momenta of the partons as relevant degrees of freedom. The inte- 
gration over transverse momenta is done up to a cutoff, /ijR, all momenta below which 
are included in a so-called hadron distribution amplitude 0, whereas larger transverse 
momenta are calculated in perturbation theory. The correlation function is hence decom- 
posed, or factorized, in perturbative contributions T and nonperturbative contributions 
0, which both depend on the longitudinal parton momenta and the factorization scale 
/iiR. If the TT is an effective quark-antiquark bound state, as is the case to leading order 
in the light-cone expansion, we can write the corresponding longitudinal momenta as up 
and (1 — u)p, p being the momentum of the vr. The schematic relation (1) can then be 
written in more explicit form as 



n+(g2,p|) = V / duT('^)(n,g2,p|,/iiR)0(")(u,/iiR). 
„ Jo 



(7) 



As n_|_ itself is independent of the arbitrary scale /ijR, the scale-dependence of T^"^ and 
(^(") must cancel each other. ^ If 0^"^ describes the meson in a two-parton state, it is called 

"^If there are more than one contributions of a given twist, they wiU mix under a change of the 
factorization scale /^ir and it is only in the sum of all such contributions that the residual /iip dependence 
cancels. 



a two-particle distribution amplitude (DA), if it describes a three-parton, i.e. quark- 
ant iquark-gluon state, it is called three-particle DA. In the latter case the integration 
over u gets replaced by an integration over two independent momentum fractions, say 
«! and a2. Eq. (7) is called a "collinear" factorization formula, as the momenta of the 
partons in the vr are collinear with the vr's momentum, and its validity actually has to be 
verified. We will come back to that issue in the next section. 

Let us now define the distribution amplitudes to be used in this paper. Again we only 
quote formulas for the vr meson, those for the K and rj are analogous. All definitions 
and formulas are well-known and can be found in Ref. [16]. In general, the distribution 
amplitudes we are interested in are related to nonlocal matrix elements of type 

{0\u{x)r[x, -x]d{-x)\7r{p)) or {0\u{x)[x,vx]rG''^^{vx)X'' /2[vx, -x]d{-x)\'iT{p)) . 

X is light-like or close to light-like and the light-cone expansion is an expansion in x^; v 
is a number between and 1 and F a combination of Dirac matrices. The expressions 
[x, —x] etc. denote Wilson lines that are needed to render the matrix elements, and hence 
the DAs, gauge- invariant. One usually works in the convenient Fock-Schwinger gauge 
x'^A'f^{x)X°'/2 = 0, where all Wilson lines are just 1; we will suppress them from now 
on. The DAs are ordered by twist, i.e. the difference between spin and dimension of the 
corresponding operators. We will include DAs of twist-2 (the leading twist), 3 and 4. The 
leading-twist DA ^^r is defined as 



(0|M(x)7/,75c/(-x)|7r (p)) =ifnP^ due'^P-'' 



(m) + - mlx^A{u) 



2 rl 

mi '^ 



+ if.— x^ due<^-^gM + 0{x,x^) (8) 

Vx Jo 

with ( = 2u — 1 and p"^ = 0. The above matrix element also contains two twist-4 DAs, g^^ 
and A. The variable u can be interpreted as the momentum fraction carried by the quark 
(as opposed to the antiquark) in the meson. 

There are two two-particle twist-3 DAs, (pp and (p^, which are defined as 

{0\u{x)t^,d{-x)\nip)) = fil [ due'^P-^M^) , (9) 

{0\u{x)iaf,^'y5d{-x)\7r{p)) = --jjl{l - pl){pf,Xu-Xf,p^) / rfMe"^^'^0<^(M) , (10) 

^ Jo 

where /i^ = Uml/{mu + m^) and pi = {rriu + rridY/ml. 

The precise definitions of three-particle DAs are a bit cumbersome and given in App. B. 
The salient feature is that there is one three-particle DA of twist-3 and four of twist-4. 

Although we have introduced not less than 10 different DAs, which are all nonper- 
turbative quantities, it may seem, at first glance, that light-cone sum rules do not retain 
much predictive power. Fortunately, however, it turns out that the DAs are highly con- 
strained functions which can be analysed in the framework of conformal expansion, a 
topic being discussed in App. B. The main result is that, to next-to-leading order in 
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Figure 1: Perturbative contributions to the correlation function H. The external quarks 
are on-shell with momenta up and (1 — u)p, respectively. 

conformal expansion, which is sufficient for the accuracy we are aiming at, all 10 DAs can 
be expressed in terms of 7 independent hadronic parameters. 

This completes the definitions necessary for the calculation of formfactors. 



3 The Sum Rules 

The diagrams to be calculated to 0{as) for two-particle DAs are shown in Fig. 1. The 
quark (antiquark) is collinear with the light meson and carries momentum up ((1 — u)p). 
Quarks are projected onto the corresponding distribution amplitudes using the complete- 
ness relation 



Uadb 



- {l)ba{ud) - - {il5)ba{ui'y5d) + - {lt,)ba{u'y''d) - - {^i ^-ir,)ba{uYl^>d) 

+ ^ {cr^^)ba{ua^"'d) . 



The diagrams are calculated in momentum space. The terms in x^ in the contribution of 
0O-, Eq. (10), are rewritten in terms of derivatives 



d 



X. 



d{up)^ 



In the previous section we mentioned that the fact that 11 can be written in factorized 
form can not be taken for granted, but requires proof. We do not attempt to give a 
proof to all orders in a^, although that should be possible using the techniques of SCET, 
but restrict ourselves to 0{as) in twist-2, to all orders in the conformal expansion, and 
to 0{as) and leading order in the conformal expansion for twist-3. The proof essentially 
relies on the cancellation of singularities, of which there are several possible types: infrared 
and ultraviolet singularities arising from loop calculations and so-called soft singularities 
which occur when the integral over u in Eq. (7) does diverge at the endpoints. The latter 
divergences have actually posed a severe problem in early attempts to treat /^ in QCD 
factorization: in Ref. [17] only the hard gluon exchange was included, which yields a 
logarithmic divergence for the parton configuration where the u quark emerging from the 



weak decay carries essentially all pion momentum. As we understand now, this divergence 
disappears when contributions from the Feynman-mechanism are added. In our case, it 
turns out that all T are regular at the endpoints m = 0, 1, so there are no soft divergences, 
independent of the end-point behavior of the distribution amplitudes. As for infrared and 
ultraviolet singularities, they can be treated in dimensional regularisation. Using the 
lowest-order expression of the Brodsky-Lepage evolution kernel for 0^ derived in [12], we 
have followed the strategy outlined in [18] to check that the infrared divergences precisely 
cancel those contained in the bare DA (j)^^^^. As for twist-3, the evolution kernel is not 
known, so we have only checked the cancellation of infrared divergences of the lowest 
order term in the conformal expansion, whose divergent behavior is well known - in fact, 
only the one-loop renormalisation of the quark condensate is needed. The ultraviolet 
divergences cancel for /_|_ and /o, which are physical formf actors and hence do not depend 
on the ultraviolet renormalisation scale; for /y, we reproduce the well-known one-loop 
anomalous dimension. 

We then have used the explicit expressions for the twist-2 and 3 two-particle DAs 
given in App. B to perform the integration over u analytically. Actually it is not the 
correlation function 11 itself that is needed, but its imaginary part, see below. 11 has 
a cut in p\ starting at rra^ and taking the imaginary part after integration over u is 
straightforward. The strategy outlined here is different from the procedure we followed 
in our previous papers [3, 4], where we took the imaginary part before integrating over u. 
This latter procedure resulted in expressions with a very complicated analytical structure 
which made it impossible to give explicit formulas for the imaginary parts. With our new 
procedure we obtain lenghty, but not very complicated expressions; the complete set of 
spectral densities p = (Imn)/7r for the sum rule for the formfactor /+ is given in App. C 

Armed with the spectral densities, we can derive the LCSR for e.g. the formfactor /_|_. 
The basic quantity is 11+ , which is calculated in two ways. In light-cone expansion, it can 
be written in dispersive representation as 



^'f{pW)= ds^-^--^ (11) 

Jml ^ Pb 

with the explicit expression for the spectral density p+^(s) given in App. C This expres- 
sion has to be compared to the physical correlation function, which also features a cut in 
p\, starting at m\: 

nf^^(p|,g^)= / dsP^^^^- (12) 

Jm^ ^ Pb 

the spectral density is given by hadronic contributions and reads 

pf'is, q') = fBmlf\{q')5{s - m,|) + pf ^~ ---(,, q^). (13) 

Here /^ is the B meson decay constant defined as 

(0|g7M75^|-B) = «/i?p^ or (mb + mg)(0|gz756|5) = m|/B. (14) 

To obtain a light-cone sum rule for /+, one equates the two expressions for n_|_ and uses 
quark-hadron duality to approximate 

^higher-mass states^^^ ^2) ^ ^LC(^^ ^2)q^^ _ ^^^^ (^5) 



where sq, the so-called continuum threshold is a parameter to be determined within the 
sum rule approach itself. In principle one could now write a sum rule 

and determine /+ from it. However, in order to suppress the impact of the approximation 
(15), one subjects both sides of the equation to a Borel transformation 

B — ^ = -^ exp{-s/M^) 



s — p% s ^ Vb ^^ 

which ensures that contributions from higher-mass states be sufficiently suppressed and 
improves the convergence of the OPE. We then obtain 



-ml/Af^^l 



(•so 

fs UiQ')= / dse-^/'^'p'fis^q'). (16) 

J ml 



This is the final sum rule for /+; expressions for the other formfactors are obtained 
analogously. The task now is to find sets of parameters M^ (the Borel parameter) and Sq 
(the continuum threshold) such that the resulting formfactor does not depend too much 
on the precise values of these parameters; in addition the continuum contribution, that is 
the part of the dispersive integral from sq to oo that has been subtracted from both sides 
of (16), should not be too large, say less than 30% of the total dispersive integral. 

4 Numerics 

In this section we obtain numerical results from the sum rules (16). The section is organ- 
ised as follows: in Sec. 4.1 we explain how we determine the sum rule specific parameters, 
i.e. the Borel parameter M^ and the continuum threshold sq. We also determine fs, 
which is a necessary ingredient in (16). In Sec. 4.2 we explain in more detail how we fix 
the hadronic input parameters, in particular the Gegenbauer moments ai^2,4 that describe 
the final state mesons. In Sec. 4.3 we calculate the formfactors at g^ = and discuss their 
uncertainties. In Sec. 4.4 we present the formfactors for central input-values of the pa- 
rameters and provide a simple parametrization valid in the full kinematical regime of q^. 
The results for g^ = are collected in Tab. 2 and Eq. (27), central results for arbitrary g^ 
in Tab. 3. More detailed results that allow one to determine the formfactors for arbitrary 
values of m^ and the Gegenbauer moments ai^2,4 are collected in App. A. 

4.1 Fixing the Borel Parameter and the Continuum Threshold 

We illustrate our procedure to determine M^ and sq with the comparatively simple ex- 
ample of Jb, the B decay constant defined in (14). This example is actually of immediate 
practical use, as Jb enters our determination of the formfactors from Eq. (16). Since it is 
not known from experiment, its value has to be taken from theoretical calculations ~ which 
basically means either lattice determinations [19] or (local) QCD sum rules [20, 21]. To 
ensure consistency of our calculations, we use the values of /b as determined from QCD 
sum rules to 0{as) accuracy [20]. The reason for this choice is twofold: firstly, it is 

8 



well-known that the use of J'b from sum rules reduces the dependence of the formfactors 
on input-parameters, in particular nib [1]; secondly, O(a^) corrections to /^ turn out 
to be rather large [21], which was anticipated in the second reference in [20], where it 
was argued that these corrections are dominated by Coulombic corrections. Precisely the 
same corrections also enter the light-cone expansion of the correlation function 11, but 
will largely cancel in the ratio /+ ~ 11/ Jb- In conclusion, we expect a cancellation of 
both large radiative corrections and parameter dependence in the formfactors when /^ 
is replaced by its sum rule to 0{as) accuracy; we do not expect the resulting numerical 
values of j'b to be "good" predictions for that quantity. 
The sum rule for Jb reads [20]^ 



flmle-^ = / dsp^''^\s)e~Jl^ + Cg,{qq) + C,G,{q<^9Gq) ^ / ds p'"\s)e-^ .{17) 



'0 



The C are the Wilson coefficients multiplying the condensates, for which we use the 
following numerical values at yU = 1 GeV: 

{qq) = -(0.24 ± 0.01)^ GeV^ and {qagGq) = 0.8 GeV'''{qq). (18) 

The condensates (and as) are actually evaluated at the scale M^. The criteria for deter- 
mining M^ and sq are often not stated very precisely. In the present context, with many 
different formfactors to calculate, which entails the need for a well-defined procedure to 
determine the input-parameters for each of them, we decide to opt for a precisely defined 
method to fix the pair (M^, sq) and impose the following criteria on the sum rule for Jb 
(and, later on, the formfactors): 

• the derivative of the logarithm of Eq. (17) with respect to 1/M^ gives a sum rule 
for the B meson mass niB'- 

so /"SO 



m%= / dssp*°*(s)/ / rfsp*°*(s). 

We require this sum rule to be fullfiUed to high accuracy ~ 0.1%. 

• the sum rule for Jb is required to exhibit an extremum for a given pair (M^, sq). 

These criteria define a set of parameters for each value of ?nfe, which are collected in Tab. 1, 
together with the resulting Jb- For all these parameter sets the continuum contribution 
(i.e. the integral J^) is between 25% and 30% of the B contribution and hence well under 
control. 

For the formfactors /^, f^ and /^ we follow the same procedure which results in 
different values of M^ and Sq for formfactors and Jb- For K and t] we use the same 
values for the Borel parameter and the continuum threshold. From the explicit formulas 
of the tree-level sum rules for the formfactors quoted in e.g. the 3rd reference in [1], one 
finds that the effective Borel parameter is uM^q rather than M^q.^ In order to keep this 

■'The contribution of the gluon condensate is not sizable and we therefore neglect it. 
^We denote the Borel parameter of the LCSR (16) by Af^c and the Borel parameter of the SR (17) 
by M2. 





nib 


So 


Af' 


Ib 


5g ~ ig 




s'o 


cf 


set 1 


4.85 


33.8 


3.8 


0.150 


33.3 


2.00 


33.6 


2.4 


set 2 


4.80 


34.2 


4.1 


0.162 


33.9 


2.25 


34.3 


2.5 


set 3 


4.75 


34.6 


4.4 


0.174 


34.5 


2.50 


35.1 


2.6 


set 4 


4.60 


35.7 


5.1 


0.210 


36.8 


3.00 


37.8 


2.9 



Table 1: Parameter sets for fs and /(O); we use the same values of Cc and sq for vr, K 
and 7]. rrih and /b are given in GeV, Sq and M^ in GeV^. 



product constant, we rescale the Borel parameter by (u) by 



{n){q')^ 



"i(,-(l-")9 



du u (t).„{u)e ^^aP / / du (j).„{u)e ^^aj^ 

no '''"0 



m2„(^_„)^2 



uo 



resulting in the approximate values (M)(OGeV^) = 0.86 and (M)(14GeV^) 
metrising the relation between the Borel parameters by 

we obtain the values and continuum thresholds given in Tab. 1. 



2 2 

mi — q 

2~' 

so-r 



0.77. Para- 



(19) 



4.2 Hadronic Input Parameters 

The hadronic parameters needed are, for each meson, 7 parameters characterising the 
twist-2, 3 and 4 distribution amplitudes to NLO in the conformal expansion, cf. App. B, 
the decay constants of the tt, K and rj and B, the factorization scale /iiR, the h quark 
mass rrib and the strong coupling a^. As for the latter, we fix as{mz) = 0.118 and 
use NLO evolution down to the required scale. The quark mass parameter entering our 
formulas is the one-loop pole mass tti^ for which we use nib = (4.80 ± 0.05) GeV (cf. 
Table 6 in the recent review [6] and references therein). We also include results for 
TTib = 4.6 GeV. The infrared factorization scale separating contributions to be included 
in DAs and perturbatively calculable terms is chosen to be /ijj^ = m\ — ml, which also 
sets the scale of a^; we will discuss the residual scale-dependence of our results below. 
The decay constants for the vr and K are very well known experimentally; for the rj the 
situation is complicated due to 77-77' mixing. We use the following values: 



/^ = 131 MeV, Jk = 160 MeV, /^ = 130 MeV. 



(20) 



fs has been discussed in the previous subsection. 

As for the meson DAs, we quote the preferred values for the twist-3 and 4 parameters 
in Tab. G; the form factors are not too sensitive to their precise values. The situation is 
different, however, for the Gegenbauer moments ai^2,4{f^) parametrizing the twist-2 DAs 
'P-K,K,ri-, and so we shall discuss in a bit more detail what is presently known about these 
parameters. 

Both theoretical calculations and experimental determinations focus mainly on the 
Tx DA (for which all odd Gegenbauer moments vanish due to G-parity; in particular 
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aj = 0). The probably earliest calculation of the lowest Gegenbauer moment a2 was done 
by Chernyak and Zhitnitsky (CZ), yielding [22] 



,cz 



(0.5 GeV) = 2/3. 



This result was obtained from local QCD sum rules, where a„ is extracted from the 
correlation function of the (local) interpolating field u'^ul^iD ■x)"'d, where x defines the 
light-cone, a;^ = 0, and the usual interpolating current for the tt, u'^^'^^d. The price to 
pay for the expansion of an intrinsically nonlocal quantity like (f).,^ in contributions of local 
operators is an increasing sensitivity to nonperturbative effects, i.e. the precise values of 
the condensates. As the coefficients of the condensates in the sum rule for a^ increase with 
powers of n and, for sufficiently large n, dominate over the perturbative contributions, it 
is clear that this method is inappropriate for calculating high moments, but one might 
expect it to be reliable at least for the lowest moment with n = 2. 

The DA obtained by CZ has the remarkable feature that 0^(1/2, 0.5 GeV) = 0, which 
is of course an artifact of neglecting all contributions from a„>4. It was subsequently 
shown by Braun and Filyanov (BF) [23] that both the pion-nucleon-nucleon coupling 
QttNN and its mesonic equivalent gpujn, when calculated from LCSRs, require a value of 
0^(1/2) significantly different from (albeit at a slightly different scale): 

3 9 45 

0^(1/2, 1 GeV) = 1.2 ± 0.3 = - - - a2(l GeV) + — 04(1 GeV) + . . . , (21) 

2 4 lb 

where the dots denote neglected terms in a.„>6. The large error is due to a large sensitivity 
of this result to twist-4 corrections to the sum rules. BF also redetermined 02, using the 
same procedure as CZ, and combining their result, which is consistent with a^^, with the 
above constraint from 07^(1/2), they obtained 

a^^{l GeV) = 0.44, a?^(l GeV) = 0.25. 

An alternative calculation aims to cure the problem of increasing condensate contri- 
butions by resumming them into nonlocal condensates [24]. The Gegenbauer moments in 
this approach are mostly sensitive to the ratio 

A2 = {qagGq)/{2{qq)) = (0.4 ± 0.1) GeV^ (/i = 1 GeV) 

and have moderate to small values. The most recent paper on that topic, Ref. [25], quotes 

02(1. 16 GeV) = 0.19, 04(1.16 GeV) = -0.13, ae.s.io ~ 10"^ (22) 

There are not too many lattice calculations of moments of the vr DA. The fairly old 
values quoted in [26] for the 2nd moment suffer from large uncertainties. This quantity has 
been investigated again recently [27], but the results, obtained in quenched approximation, 
are still preliminary. 

Alternative determinations of Gegenbauer moments rely on the analysis of experimen- 
tal data, in particular the pion-photon transition formfactor 7 + 7* -^ vr, measured at 
CLEO and Cello, and the electromagnetic formfactor of the pion. The results of these 
analyses are typically either determinations of 02 (setting a„>4 to 0) or constraints on 
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a linear combination of 02 and a^ (setting an>6 to 0)7 These determinations are lim- 
ited by mainly two problems: large experimental errors and the contamination by poorly 
known twist-4 and higher effects, which are usually estimated from QCD sum rules. As 
for the pion-photon transition formfactor, which has been measured by CLEO and Cello, 
the technique used to extract 02 and 04 has been pioneered by Khodjamirian [29], re- 
fined by Schmedding and Yakovlev [30], with subsequent further refinements by Bakulev, 
Mikhailov and Stefanis [31]. The upshot is that for not too small Q^ the pion-photon 
transition is mostly sensitive to a like-sign combination of 02 and 04. Summarizing the 
analyses of this process, we conclude from Tab. I in [25] that 

02(1 GeV) + 04(1 GeV) = 0.1 ±0.1 (23) 

is a fair reflection of the current state of knowledge of 02,4 from that process. 

As for the pion electromagnetic formfactor, the authors of Ref. [32] unfortunately only 
obtain a value for 02 and set 04 to zero. A very recent analysis of that formfactor, Ref. [25] , 
concludes that calculations using the nonlocal-condensate model are in good agreement 
with data. 

So what then do we actually know about 02 and 04? It seems to us that, taking ev- 
erything together, and with due consideration of the respective strengths and weaknesses 
of different approaches, the most reliable constraints for these quantities are (21) and 
(23). These two constraints contain opposite-sign combinations of 02 and 04 and hence 
are about equally sensitive to both parameters. The resulting allowed area for 02 and 04 
is shown in Fig. 2; its center is at 

02(1 GeV) = 0.115, 04(1 GeV) = -0.015, 
02(2.2 GeV) = 0.080, 04(2.2 GeV) = -0.0089. 

These are the central values we will use in our calculation of formfactors. The figure also 
shows that the remaining uncertainties are still considerable. Anticipating a future better 
determination of these parameters, from lattice or else, we will present our final results 
in such a way as to facilitate the inclusion of any shift in these values. Since much less 
is known about the Gegenbauer moments of the other pseudoscalar mesons K and rj, we 
resort to SU(3) symmetry and use the same Gegenbauer moments. 

Eq. (24) and Fig. 2 confirm the findings of previous analyses that the CZ DA is strongly 
disfavored; the same applies to the values obtained by BF and to the local QCD sum rule 
for a2, which favors a large positive 02 ~ 0.4. One explanation for the failure of the 
corresponding QCD sum rule could be that already the case n = 2 may be too "nonlocal" 
for sum rules to work. Another one could be that the treatment of ai and other resonances 
contributing to that sum rule may be insufficient. We leave a further discussion of that 
question to future work. The result from sum rules with nonlocal condensates [24, 25, 31], 
shown as black square in Fig. 2, is also outside the favored area in Fig. 2, which is mainly 
due to the large value of |a4|. It would definitely be very interesting to see all these results 
and constraints on 02 4 be supplemented by lattice determinations. 



^In principle it is possible to determine 02, 04 and even higher moments separately from the Q^ 
dependence of their respective contributions. However, such an analysis requires accurate measurements 
of the formfactors over a large enough range of Q^, which are presently not available. Sec also Ref. [28], 
in particular Fig. 4. 
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Figure 2: 02(1 GeV) and 04(1 GeV) as determined from the constraints (21) and (23). 
Solid line: central values, dashed lines: uncertainties. The black square labeled BZ 
denotes the central values used in this paper, Eq. (24), BMS the prediction of the nonlocal 
condensate model, Eq. (22), rescaled to /z = 1 GeV, and BF is the central value obtained 
in Ref. [16]. 

The only parameter left to discuss is a\ for the K meson (by which we understand 
an sq bound state), which is a G-parity breaking parameter. Here the situation is even 
worse than for 02,4, as neither size nor even sign of that quantity are reliably known. 
The facts at hand are the following: the intuitive expectation is that a\ (i.e. the moment 
with a weight-function proportional to 2m — 1) should be positive, as the DA is expected 
to be slightly tilted towards larger values of u which is the momentum fraction carried 
by the (heavier) s quark in the meson - the heavier the quark, the more the DA is 
expected to peak at large u, the extreme case being a hq bound state whose DA should 
be close to 5(1 — u). The (tree-level) QCD sum rule calculation in [22] seemed to confirm 
intuition, but was challenged, when Ref. [33] found a sign-mistake in that calculation and, 
including two- loop radiative corrections, obtained a negative sign for a^ . For this paper, 
we first decided to stick to that result and use the central value af(l GeV) = —0.18. It 
turned out, however, that this value tends to produce formfactors with an unfavorable 
g^-dependence.^ We therefore decided to revert to the original result by CZ [22] and use 



K 



(1 GeV) = 0.17 



K 



(2.2 GeV) = 0.135. 



(25) 



The conclusion from that inconclusive situation can only be that a second opinion has to 
be sought, and we urge our colleagues from the lattice community to take up the challenge 
and provide the first-ever lattice determination of af^. For the time being, we will present 
our results in a way that makes it possible to obtain the formfactors also for different 
values of a^ . 



4.3 



Results for o^ = 







Let us first discuss the sum rule results for q^ = 0. They are collected in Tab. 2, for all 
4 parameter sets from Tab. 1.^ Including the uncertainty of mi,, rub = (4.80 ± 0.05) GeV, 
the final central values and uncertainties of the formfactors are given in Eq. (27). 

^That is: formfactors not very compatible with the paramctrization discussed in Sec. 4.4, which is 
based on generic analytic properties of the formfactors. 
Vo(0) is not included as /o(0) = /+(0). 
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Figure 3: Central values of the formfactors /(O) and uncertainties A. Numbers from 
Tab. 2. 
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Figure 4: Dependence of /.^(O) on 02 (/xm) and a4(/iiR), for parameter set 2. The lines 
are lines of constant /+(0). The dot labeled BZ denotes our preferred values of 02,4, BMS 
the values from the nonlocal condensate model and BF from the sum rule calculations of 
Ref. [16]. 
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Figure 5: (a) Dependence of /:f^(0) on the Gegenbauer moment ai(/iiR). (b) f^iq"^) as 
function of g^ for different values of ai: solid line: af^(lGeV) = 0.17, short dashes: 
af (IGeV) = 0, long dashes: af^(lGeV) = —0.18. Input parameters: set 2. 
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set 1 


set 2 


set 3 


set 4 


Aas 


^a2,a4 


A 


A„, 


fm 


0.250 


0.258 


0.263 


0.274 


0.023 


0.019 


0.030 


— 


fm 


0.244 


0.253 


0.260 


0.273 


0.013 


0.022 


0.026 


— 


/f(0) 


0.324 


0.331 


0.335 


0.339 


0.033 


0.023 


0.040 


0.256a, 


/#(0) 


0.347 


0.358 


0.367 


0.381 


0.022 


0.027 


0.035 


0.316a, 


fm 


0.269 


0.275 


0.278 


0.286 


0.029 


0.019 


0.035 


— 


fm 


0.277 


0.285 


0.292 


0.305 


0.018 


0.022 


0.028 


— 



Table 2: Final central values of the formfactors at g^ = for the parameter sets of Tab. 1. 

/o(0) = /+(0). The errors A„,, A 

A = (A^ 

on the sign of a^ and can become negative. 



as, '-'02,04 and Aa^ are described in the text. A is defined as 
A^2 aiY^'^ ^^d 6a, as 6a, = ai(lGeV) — 0.17. Note that 6a, carries information 



The formfactors are calculated from Eq. (16) using the parameter sets given in Tab. 1 
and the hadronic input parameters given in Eqs. (24) and (25) and Tab. G. The depen- 
dence of the formfactors on mi,, i.e. the set, is shown in Fig. 3. It is evident that the 
residual dependence of /(O) on mb is much smaller than the one of Jb in Tab. 2, which 
confirms our expectation that the calculation of Jb from a sum rule reduces the parameter 
dependence of the formfactors. /+(0) depends sensitively on 02 and 04 as illustrated in 
Fig. 4. The formfactors show moderate SU(3) breaking between vr and 77, which is due to 
terms in the LCSRs proportional to the meson mass. For K, the situation is different, and 
we observe a strong enhancement of the formfactor due to the combination of two effects: 
the fact that fx is larger than f^^ and the positive contribution of the Gegenbauer moment 
fli to the formfactor. As discussed in the previous subsection, the numerical value of Oi, 
and even its sign, is not precisely known. Fig. 5(a) illustrates the dependence of /^(O) 
on Oi, which is quite strong. Fig. 5(b) shows the dependence of f^{q^) on q^ for different 
values of ai. It is evident that ai mainly determines the normalisation of the formfactor, 
but has only minor impact on its shape. The uncertainty of /_f^(0) induced by ai will 
be discussed below. The dependence of /+(0) on the sum rule parameters M^ and sq is 
illustrated in Fig. 6 and is very mild, thanks to the optimized criteria for choosing M^ and 
So outlined in Sec. 4.1. The behavior of the other formfactors is very similar. In Fig. 7 we 
show the variation of /+(0) with a change of the factorization scale /xir in the large range 
1 GeV < /iiR < mi). The curve is remarkably fiat which can be understood from the fact 
that radiative corrections cancel to a certain extent between 11+ and fs and that large 
logarithms of type himi,/ ^m_ occur only at subleading order in the conformal expansion of 
the DAs, which is numerically suppressed with respect to the leading (/iiR-independent) 
term, and at subleading twist, which is also suppressed. 

Let us now turn to the uncertainties of the formfactors induced by a variation of 
the input parameters. It is convenient to split the formfactors into contributions from 
different Gegenbauer moments. 



f{q') = r(g') +«ir (?') + {a2r (g^) +a4r(g2)} 



(26) 



where /"* contains the contributions to the formfactors from the asymptotic DA and also 
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Figure 6: Dependence of /+(0) on (a) the Borel parameter M^ and (b) the continuum 
threshold Sq. Input parameters: set 2 in Tab. 1. 
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Figure 7: Dependence of /+(0) on the factorization scale /ijR. Same input parameters as 
in Fig. 6. 

all higher-twist effects from three-particle quark-quark-gluon matrix elements. Explicit 
expressions for the functions y'^«'"i'"2,a4 f^g^^ ]^g obtained from Tab. C in App. A; in par- 
ticular /"*(0) is just given by the parameters a in that table. We calculate separately the 
uncertainties Aas,ai of the first and second term and the combined uncertainty A^j^a^ of 
the term in curly brackets. We start with A^^. To estimate its value we vary the following 
quantities: 

• the threshold sq by ±0.5 GeV^; 

• the Borel parameter M^ in Eq. (19) by ±1.2 GeV^; 

• the infrared factorization scale fif^ = m\ — ml by ±2 GeV^; 

• the quark condensate and the mixed condensate as indicated in Eq. (18); 

• the twist-3 matrix-element 773 by ±50%. 

nib is kept fixed and we calculate the uncertainty separately for each parameter set; for 
a given formfactor, Aas is then the largest uncertainty of the 4 sets. The errors are 
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correlated and we therefore scan the five-parameter space for the largest deviations from 
the central values. The resulting A^^ are given in Tab. 2. 

The uncertainty of f^{0) induced by ai is dominated by ai itself, so we do not attempt 
to determine the uncertainty of /"^ arising from varying M^, sq etc., but just take the 
maximum value of /"^(O) = a from Tab. C in App. A and multiply it by 6i = ai(l GeV) — 
0.17 and the leading-order scaling factor from IGeV to /ijR, which gives the entry labeled 
Aq^ in Tab. 2. 

As the allowed input values of 02 and 04 are correlated and given by the rhomboid 
shown in Fig. 2, we only determine the combined uncertainty Aa2,a4, arising from the 
corresponding variation of the Gegenbauer moments, separately for each parameter set. 
The resulting uncertainties depend strongly on the precise values of M^ and Sq, so for a 
conservative estimate of the uncertainty we scan the full 7-parameter space in 02, 04, M^ 
etc. and quote the largest deviation from the central value as uncertainty, which yields 
the Aa2,a4 quoted in Tab. 2. Taking everything together, and including the variation of 
rrif, = (4.80 ± 0.05) GeV in the error estimate, adding errors in quadrature, we find ((5a^ is 
defined in Tab. 2): 



fm 


= 0.258 ±0.031, 


fm 


= 0.253 ±0.028, 


/^(O) 


= 0.331 ± 0.041 + 0.25(5a„ 


/f(0) 


= 0.358 ± 0.037 ±0.31<5„„ 


/:(o) 


= 0.275 ±0.036, 


im 


= 0.285 ±0.029. 



(27) 

These are our final results for the formfactors at g^ = 0. For f"^'^ the total theoretical 
uncertainty is 10% to 13%, for f^ it is 12%, plus the uncertainty in ai, which hopefully 
will be clarified through an independent calculation in the not too far future. These 
uncertainties include a variation of both the external input parameters and the sum rule 
specific parameters, but they do not include an additional "systematic" uncertainty of 
the sum rule method itself. To a certain extent, this intrinsic sum rule uncertainty is 
included by the variation of the sum rule specific parameters M^ and sq, which sets the 
minimum uncertainty of the result: all external hadronic parameters fixed, this variation 
induces a ~ 7% uncertainty of /+(0) quoted in Eq. (27). Realistically, one may hope to 
reduce the ~ 12% uncertainty quoted to ~ 10% by reducing the errors on the Gegenbauer 
moments 02,4 by a factor of 2. Further improvement will then have to come from a better 
control over higher-twist matrix elements, dominated by the quark condensate and the 
quark-quark-gluon matrix element 773 discussed in App. B. 

4.4 Results for q^ 7^ 0, Fits and Extrapolations 

In this subsection we calculate the g^-dependence of the formfactors for central values of 
the input parameters and cast them into a three-parameter parametrization that is valid 
for all q^. The results are given in Tab. 3 which is to be used together with Eq. (30). 
The fit parameters for other sets of input parameters are given in App. A. We refrain 
from a complete analysis of the uncertainty of the g^-dependence of the formfactors, but 
just mention that it is likely to be smaller than that at g^ = 0, which is indicated by a 
decrease of the spread between the formfactors calculated from the different parameter 
sets in Tab. 1, cf. Fig. 9. 
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The validity of the LCSR approach is restricted to the kinematical regime of large 
meson energies, Ep ^ Aqcd, which via the relation 



q = nig — 2mBEp 

implies a restriction to small and moderate g^; specifically, we evaluate the sum rules 
only for < g^ < 14GeV^. The resulting formfactors are plotted in Fig. 8, using the 
parameter set 2 in Tab. 1 and the hadronic input parameters given in Eqs. (24) and (25) 
and Tab. G. As expected from LEET [34], /+ and /t nearly coincide. Although this 
agreement is expected to be best for small g^, i.e. large energies of the light meson, it is 
seen to hold for all q^. From the LCSR point of view, this agreement is due to the fact 
that the leading twist-2 contributions to the corresponding correlation functions coincide 
at tree-level. The figure also shows that the g^-dependence of /o is weaker than that of the 
other formfactors. This is can be understood from the fact that, if /+ is represented as a 
dispersion relation over hadronic states, these states have quantum numbers J^ = \~ and 
hence zero orbital angular momentum, whereas for /q the quantum number is J^ = O"*" 
and thus the coupling of these states or, in the language of potential-models, their wave- 
function at the origin, is suppressed as it corresponds to a state with orbital angular 
momentum L = 1. Fig. 8 also shows sizable SU(3) breaking for the K, but a moderate 
one for r], which is due to the same effects discussed for the formfactors at g^ = 0. In Fig. 9 
we show /^(g^) as function of g^, calculated for sets 1, 3 and 4 and normalized to set 2. It 
is evident that the uncertainties induced by rrif,, which amount to 6% at g^ = 0, become 
less important for larger g^, so that for instance the branching ratio of the semileptonic 
decay B —>■ nev will be less dependent on the precise value of rrt}, than /+(0). 

One of the main goals of this paper is to give simple expressions for the formfactors 
in the full physical regime < g^ < (m^ — mpY ~ 23 GeV^. We thus have to find a 
parametrization that 

• reproduces the data below 14 GeV^ with good accuracy; 

• provides an extrapolation to g^ > 14 GeV^ that is consistent with the expected 
analytical properties of the formfactors and reproduces the lowest-lying resonance 
(pole) with J^ = 1~ for /+ and /t-^° 

It is actually not very difficult to find good fits: the parametrization 

fi<l') = . , ,..,i^°l, ..,.,2 (28) 



1 — apq^ /m\ — hp^q^ / 



advocated in previous works, e.g. [4], is one example for an excellent fit to the results 
of the sum rules for g^ < 14GeV . In the present context, however, it turns out to be 
unsuitable as it produces, for /^, a pole at g^ ~ 23 GeV^, which is below the physical 
pole at g^ = m^. = (5.32 GeV)^. In our previous paper [4] we argued that the above 
parametrization should be matched to a simple pole- dominance formula /+ ~ l/(m^* —(f) 
for g^ above a certain threshold gg ~ 15 GeV^, defined as the value of g^ that would allow a 

lOpor /o; the lowest pole with quantum numbers O'^ lies above the two-particle threshold starting at 
{niB + mpY' and hence is not expected to feature prominently. 
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Figure 8: /+ (solid lines), /o (short dashes) and fx (long dashes) as functions of g^ for vr, 
i^ and f]. The renormalisation scale of /t is chosen to be rrif,. Input parameters: set 2 in 
Tab. 1. 



smooth transition^^ from one parametrization to the other. This procedure unfortunately 
does not work for our new formfactors, as the optimum Qq turns out to be far outside the 
physical regime. We therefore decide to follow, as far as possible, the procedure advocated 
by Becirevic and Kaidalov [35], who suggested to write the formfactor /+ as dispersion 
relation in g^ with a lowest-lying pole plus a contribution from multiparticle states, which 
in turn is to be replaced by an effective pole at higher mass: 



/+(g^ 



ri 



1 — q^ /m\ 



ds 



{niB+mp)'^ 



PJs) 
s — q^ 



ri 



q^ /m\ 



+ 



^2 



q^ /m 



(29) 



(30) 



The lowest-lying resonance in the hu channel is well known experimentally: it the B*{1~) 
vector meson with mass 5.32 GeV; this is also the mass to be used for the 77, as the 
B —>■ 7] formfactors calculated in this paper refer to a. b —>■ u transition. For the K we 
have calculated the mass of the B* resonance in the heavy-quark limit and find 



m 



Bi 



— m. 



m. 



— m 



B 



m 



K 



niB* = 5.41 GeV. 



For Eq. (30) to describe all /+ and also J't, which feature the same 1 resonance, 
in terms of three fit parameters, ri, r2 and rragt, it is crucial that the position of the 
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That is equality of both the parametrization formulas and their first derivatives in Qq . 
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Figure 9: Ratio of /^ yO^)/ f+ i'f) &s function of g^. Solid line: set 1; long dashes: 
set 3; short dashes: set 4. 
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Figure 10: Variation of the total semileptonic rate r(i? — *• i^ev) as function of g'^ax' ^^^ 
maximum q^ for which LCSR results are included in the fits. The rate is normalized to 
1 for g^a^ = 14GeV^ and fit 1. Solid line: fit 1, long dashes: fit 2, short dashes: fit 3. 
Input parameters: set 2. 

lowest pole is sufficiently below the two-particle cut starting at (m^ + mp)"^. We find 
that indeed most /^ rp formfactors, with the exception of /^ , are described very well 

by (30). For /^ , however, and all /_,_^, met gets too close to mi, so that the fit 
becomes numerically unstable. In this case, it is appropriate to expand (30) to first order 
in TTifit — rrii , which yields 



f^:?i-'' 



\q' 



n 



q^ /m\ 



+ 



r-2 



q^ /m\ 



2N2 



(31) 



with fit parameters ri and r2, and irti = mB*^B* fixed. 

For /o, one can write a decomposition similar to (29), but here the lowest-lying pole 
with quantum numbers O"*" lies either above the two-particle threshold (for vr and f]) or is 
very close to it (for K, cf. Tab. D), so that the pole is effectively hidden under the cut 
and only the dispersive term survives in (29). We again follow the suggestion of Becirevic 
and Kaidalov and replace this term by an effective pole, i.e. we set 



fo{q' 



^2 



q'^/m 



(32) 



The accuracy of the fits of the LCSR results to the above parametrizations is generally 
very high and best for sets 1 to 3 of Tab. 1 with mb = (4.80 ±0.05) GeV, with a maximum 
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ri 


r2 


{m,Y 


Ht 


n 


0.744 


-0.486 


{mlf 


40.73 


fS 





0.258 


— 


33.81 


fr 


1.387 


-1.134 


{m\f 


32.22 


n 


0.162 


0.173 


«)' 


— 


/o^ 





0.330 


— 


37.46 


f¥ 


0.161 


0.198 


{raff 


— 


n 


0.122 


0.155 


K)2 


— 


fo 





0.273 


— 


31.03 


Pt 


0.111 


0.175 


{mlf 


— 



Table 3: Fit parameters for Eq. (30) for set 2 in Tab. 1 and central values of the input 
parameters of the DAs, Eqs. (24), (25) and Tab. G. mi is the vector meson mass in the 
corresponding channel: m^''' = m^* = 5.32 GeV and mf = m^. = 5.41 GeV. The scale 
of /t is yu = 4.8 GeV. 





set 1 set 2 set 3 set 4 


fit 1 


0.97 1 1.01 1.05 


fit 2 


0.97 0.98 0.99 1.00 


fit 3 


0.95 0.98 1.00 1.04 



Table 4: Total semileptonic decay rates r{B -^ neu) normalised to 1 for set 2, fit 1, for 
different formfactor parametrizations and input parameter sets. 

1.2% deviation; set 4 fares slightly worse with an accuracy of 2% or better. The quality 
of the fits is discussed in more detail in App. A. The uncertainty introduced by fitting 
is much smaller than the actual uncertainty of the sum rule calculation, which we have 
found to be around 10% at g^ = 0, and also much smaller than the intrinsic and irreducible 
sum rule uncertainty, which we have estimated to be ~ 7%. Nevertheless it is legitimate 
to ask whether the extrapolation of the fits to g^ > 14GeV^, or the variation of the 
"cutoff" g^a^x = 14GeV^, introduce an additional uncertainty. In answering this question, 
we first would like to recall that for most applications it is actually sufficient to know the 
formfactors for g^ < 14 GeV^ only — these include in particular nonleptonic B decays 
treated in QCD factorization, and also the rare decays B -^ (tt, K, ri)£~^i~ , as the spectrum 
for invariant lepton masses above the cc threshold, i.e. g^ > 'tIj ,^ f« 10 GeV^, is dominated 
by long-distance processes unrelated to B ^ (vr, K, rj) formfactors. The only, but very 
important case where the formfactor is needed over the full range of g^ is the semileptonic 
decay B — > niu, which depends on /^ and (for decays into r) on f^. We discuss the effect 
of the extrapolation on this decay by studying three different parametrizations of /^: 

fit 1: Eq. (30), our standard parametrization; 
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fit 2: a modified version of (28), with one zero of tlie denominator fixed at ml = m\,: 



fUq' 



q'^/mf){l — q'^/m 



fit J 



fit 3: a parametrization similar to (31), but with the pole mass as fit parameter: 

^+^^^ 1 - g7< (1 - ?VO^ ■ 

We quantify the difference between these parametrizations by calculating the semileptonic 
decay rate, the integral of Eq. (4) over g^ from to (m^ — ?Ti^)^, normalizing to our 
central values, set 2 and fit 1. The results are collected in Tab. 4. It is evident that the 
dependence of the rate on the fit is rather mild, despite the double-pole of fit 3, which is 
however sufficiently far away from the endpoint of the spectrum, TTigt = (5.6 ± 0.1) GeV, 
and hence has only moderate impact on the rate. We conclude that the extrapolation 
of /^ causes an uncertainty in the total semileptonic decay rate r{B ^ ireu) which is 
considerably less than the expected intrinsic sum rule uncertainty of ~ 14%. 

We conclude the discussion of the uncertainty of the extrapolations by studying the 
effect of changing the maximum value of g^ for which the sum rules results are included 
in the fits. Our default value q"^^^ is 14 GeV^; lowering q"^^^ changes the fit parameters of 
all three parametrizations and hence the predictions for the total semileptonic decay rate. 
Fig. 10 shows the corresponding change in the rate, normalised to our central values fit 1 
and g^j^x = 14 GeV^. Again the dependence of the rate on q"^^^ is mild, which corroborates 
our conclusion that the precise shape of the formfactor is not that relevant, as long as it 
does not exhibit too strong a singularity at g^ = (5.32 GeV)^. 

There are also other tests and checks for the validity of the extrapolation of (30) to 
the full physical regime g^ < {rriB — mp)"^: firstly, the coefficient ri for /^ is related to 
the coupling gBB*n as 

n = -T, ■ 33 

At the upper end of the physical range in g^ we can expect vector-meson dominance to be 
effective and therefore the fit-parameter should be close to the above value. In fact lattice 
[36] and meson-loop calculations (cf. the first reference in [6]) yield ri ~ 0.8, but are at 
variance with a determination oi gsB'-K from LCSRs which yields ri k, 0.44 [5]. The lattice 
and meson-loop calculations are further supported by the agreement of their predictions 
for Qdd'-k with experimental measurements, whereas LCSRs again give a value that is too 
low by almost a factor of two. The author of [37] speculates that this discrepancy may be 
explained by a violation of quark-hadron duality in the LCSRs used for the determination 
of gDD*TT and gBB*-Ki which would preclude a clean determination of these couplings from 
LCSRs. Another possible solution of the problem was suggested in Ref. [38], where it was 
shown that the value of ri from LCSRs increases once a radial excitation with negative 
residue is included in the hadronic parametrization of the correlation function. ^^ If we 
interpret our fit results as determinations of gBB*TT, we get the following values of ri for 
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Note that the corresponding spectral function is not positive definite. 
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the sets 1 to 4: (0.73,0.74,0.77,0.94) (cf. Tab. A), which is in reasonable agreement with 
lattice and meson-loop calculations. 

Secondly, there is one further constraint on the formfactor /q. As first pointed out in 
Ref. [39], in the soft-pion limit p -^ and ml -^ (i.e. g^ = rra^) fQ{fn\) is related to 
the decay constants of the B and vr as 

r,{<) = ^- (34) 

In 

We can compare this relation with our parametrization by solving it for fs- For the 
four parameter sets of Tab. 1, we get from Eq. (34) /f*^ = 201 MeV, /f*^ ^ 193 MeV, 
y^et3 _ 190 MeV and f^^'^ = 207 MeV, which is in good agreement with lattice and sum 
rule calculations. 

Let us conclude with one more remark. In LEET, /_|_ and /o are related as [34]: 

/o = — /+, (35) 

TUB 

which is valid in the combined limits rriB -^ C)0 and E ^ 00. This constraint was used 
in Ref. [35] to reduce the number of fit parameters to two as necessitated by the limited 
accuracy of the lattice formfactors. We do not impose this constraint explicitly, but find 
that it is valid to 4% accuracy for our formfactors, for not too large q^. 

Summarizing, we conclude that, for all formfactors, the three-parameter formula (30) 
provides both an excellent fit to the LCSR results for q^ < 14 GeV^ and a smooth extrap- 
olation to 14GeV^ < q^ < {rriB — rrip)'^, and is consistent with all known constraints. 

5 Summary &; Conclusions 

In this paper we have given a thorough and careful examination of the predictions of QCD 
sum rules on the light-cone for the formfactors /+, /o and /^ for the decays B — ^ TT^K^rj. 
We have not discussed B — ^ 77', which is not accessible within the method due to its large 
mass. 

The main improvements of our results with respect to our previous publications [3, 4] 
are: 

• predictions for all formfactors of S -^ 7i,K,ri transitions to 0{as) accuracy for 
twist-2 and 3 two-particle contributions; 

• a well-defined and precise method for fixing sum rule specific parameters (cf. Sec. 4.1); 

• a careful assessment of uncertainties at zero momentum transfer (cf. Sec. 4.2 and 
4.3); 

• a detailed breakdown of the contributions of different Gegenbauer moments Oj to 
the formfactors (cf. App. A), which 

— renders straightforward the implementation of future updates of these param- 
eters; 
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— allows the assessment of the impact of nonasymptotic twist-2 distribution am- 
plitudes on QCD factorised nonleptonic B decays in a coherent way, to 4th 
order in the conformal expansion; 

• a parametrization of the g^-dependence of formfactors valid in the full physical 
regime of momentum transfer that reproduces all relevant analytical properties of 
the formfactors (cf. Sec. 4.4). 

Our main results for q^ = are collected in Tab. 2 and Eq. (27). They depend crucially 
on the values of the Gegenbauer moments describing the twist-2 distribution amplitudes 
oi IT, K and rj, cf. App. B. We have determined these parameters as discussed in Sec. 4.2, 
but a better determination from an independent source, e.g. lattice calculations, would be 
extremely useful. This applies in particular to the SU(3) breaking parameter af^, whose 
size and even sign is under discussion (cf. Ref. [33]). Once more precise values for these 
parameters will be available, it is straightforward to obtain the corresponding formfactors 
from the data collected in App. A. Setting aside ai, the total theoretical uncertainty of the 
formfactors at g^ = is 10% to 13%, which includes a variation of all input parameters. 
It can be further improved by reducing the uncertainties of, in particular, 02, 04, the 
quark condensate and 773, the dominant quark-quark-gluon matrix element. A reduction 
of the uncertainty of 02,4 by a factor of two will give a ~ 2% gain in accuracy, reducing 
the uncertainty of the quark condensate and 773 by the same factor will give another 2%. 
The uncertainty due to the variation of only the sum rule specific parameters is 7%, 
which cannot be reduced any further and hence sets the minimum theoretical uncertainty 
that can be achieved within this method. Comparing with the uncertainties quoted in 
our previous publications, we have achieved a reduction of the global estimate ~ 15% 
quoted in [3] and also of the 20% uncertainty for /+(0) quoted in [4]. This is partially 
due to a reduction of the uncertainties of the hadronic input parameters, in particular 
rrih, and partially due to a refinement of the assessment of sum rule specific uncertainties 
as discussed in Sec. 4.1. 

We have also calculated all formfactors for < g^ < 14 GeV^; the upper bound on g^ is 
due to the limitations of the light-cone expansion which requires the final-state meson to 
have energies E :^ Aqcd: for q^^^ = 14 GeV^ the meson energy is E = 1.3 GeV. In order 
to allow a simple implementation of our results, we have given a parametrisation that 
includes the main features of the analytical properties of the formfactors and is valid in 
the full physical regime < g^ < {mB — mpY. The corresponding results for our preferred 
set of input parameters are given in Tab. 3; a detailed breakdown of the contributions of 
different parameters to the full formfactors is given in App. A. The main features of the 
results are that the formfactors /+ and /t are nearly equal as predicted by LEET and 
that /o is very well described by a single-pole formula. The uncertainty induced by the 
extrapolation of the parametrization to larger momentum transfers is an issue only for 
the semileptonic decay B -^ neu; we have checked that the change of the total rate is at 
most 5% for three different extrapolations of the light-cone sum rule results. 

Our approach is complementary to standard lattice calculations, in the sense that 
it works best for large energies of the final state meson (i.e. small g^), whereas lattice 
calculations work best for small energies - a situation that may change in the future with 
the implementation of moving NRQCD [10]. Previously, the LCSR results for /Tq at 
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set 2, mf, 


= 4.8 GeV 




set 4, mf, = 4.6 GeV 




ri 


ml 


r2 


Ht 


A 


ri mi r2 m|t A 


n 


0.744 


{mlf - 


-0.486 


40.73 


0.3 


0.944 {mlf -0.669 34.27 0.3 


n 





— 


0.258 


33.81 


0.1 


- 0.270 33.63 1.2 


It 


1.387 


{mlf - 


-1.134 


32.22 


0.5 


use (A.2) with n = 0.152, 
r2 = 0.122, mi = m^, A = 0.4 


J+ 


0.918 


{mlf - 


-0.675 


38.20 


0.1 


0.711 {mlf -0.441 44.31 0.1 


fK,as 





— 


0.244 


30.46 


0.8 


- 0.270 31.93 0.1 


p7r,as 

It 


1.556 


{mlf - 


-1.321 


32.56 


0.2 


1.331 {mlf -1.061 33.43 0.4 



Table A: Fit parameters for the vr Eq. (A.l) for both the full formfactors and the asymp- 
totic ones, /"*, Eq. (A. 5), using the sets 2 and 4 in Tab. 1. The formfactor /q is fitted to 
the parametrization (A. 3). The mass parameters m,f are given in Tab. D. A is a measure 
of the quality of the fit and is defined in (A. 4). 

small and moderate q^ were found to nicely match the lattice results obtained for large 
q^ [42]. The situation will have to be reassessed in view of our new results and it will be 
very interesting to see if and how it will develop with further progress in both lattice and 
LCSR calculations. 
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Appendix 

A Fit Parameters and Comments 

This appendix extends the discussion of Sec. 4.4. 



Full formfactors. As discussed in Sec. 4.4, we fit the LCSR results to the following 
parametrizations: 



for fl^: 



13 



fig' 



Tl 



^2 



q'^/mf 



q'^/m? 



(A.l) 



fit 



13 



K,ri 



Apart from f^ for set 4, which shows the same behavior as / , 'j! and hence is parametrised the same 
way, i.e. according to (A.2). 
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set 2, rrib = 


= 4.8 GeV 




set 4, rrib = 


= 4.6 GeV 






ri 


r2 


mfit(mi) 


A 


n 


^2 


mfit(mi) 


A 


/^ 


0.1616 


0.1730 


i<r 


1.2 


0.1903 


0.1478 


(mf)^ 


1.0 


/o^ 





0.3302 


37.46 


1.0 





0.3338 


38.98 


1.9 


/f 


0.1614 


0.1981 


(mf)2 


0.5 


0.1851 


0.1905 


(mfr 


1.7 


fl 


0.1220 


0.1553 


K)2 


1.0 


0.1380 


0.1462 


K)2 


0.9 


r. 





0.2734 


31.03 


0.5 





0.2799 


30.46 


2.0 


Pt 


0.1108 


0.1752 


K)2 


0.5 


0.1160 


0.1841 


K)2 


1.6 


rK,as 

J + 


0.0541 


0.2166 


(mf)^ 


0.2 


0.0991 


0.2002 


(mf)2 


0.6 


pK,as 

Jo 





0.2719 


30.33 


0.7 





0.2984 


31.99 


0.5 


f-K,as 

Jt 


0.0244 


0.2590 


(mf)2 


0.8 


0.0660 


0.2621 


«? 


1.3 


rr 


0.0802 


0.1814 


K)2 


1.0 


0.1201 


0.1636 


{m^,f 


0.6 


Jo 





0.2604 


28.80 


0.5 





0.2803 


29.59 


0.8 


nrt,as 

Jt 


0.0570 


0.2115 


K)2 


0.3 


0.0914 


0.2096 


{m\f 


1.0 



Table B: Fit parameters for K and r] for Eq. (A. 2), for both the full formfactors and the 
asymptotic ones, /"*, Eq. (A. 5), using the sets 2 and 4 in Tab. 1. The formfactor /o is 
fitted to the parametrisation A. 3. The mass parameters mi are given in Tab. D. A is a 
measure of the quality of the fit and is defined in (A. 4). 

where m\ is the mass of i?*(l^), m\ = 5.32 GeV; the fit parameters are ri, r2 and 
mm; 



for Z^;;? and /? (set 4): 



IXq' 



n 



r2 



1 — q^ /m\ (1 — q^ /mi) 



2N2' 



(A.2) 



where mi is the mass of the 1 meson in the corresponding channel, cf. Tab. D; the 
fit parameters are ri and r2; 



for /o: 

fo{q') 

the fit parameters are r2 and rngt. 



r2 



q^ /mi 



2 ' 
fit 



(A.3) 



The fit parameters are collected in the upper halves of Tabs. A and B. A is a measure of 
the quality of the fit and defined as 



A = 100 max 

t 



fit) - r{t) 



fit) 



, te{0,i,...,f,14}GeV2 



(A.4) 



i.e. it gives, in per cent, the maximum deviation of the fitted formfactors from the original 
LCSR result for (f < 14 GeV^. From the A given in the table we conclude that the overall 
quality of the fits is very good and best for the pion and also that they work better for 
our preferred set 2 than for set 4. 
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set 2, m6=4.8 GeV 


set 4, mb=4.6 GeV 




a 


6x 10^ 


ex 10^ 


dx 10^ 


5 


a 


hx W 


ex 10^ 


dxW 5 


/^(oi) 


0.310 


0.930 


0.139 


-0.083 


0.3 


0.276 


0.060 


0.151 


-0.157 0.7 


/(f(«l) 


0.308 


0.106 


0.026 


-0.048 


0.2 


0.273 


-0.433 


0.0001 


-0.051 0.2 


/f(oi) 


0.381 


1.056 


0.167 


-0.108 


0.3 


0.354 


0.027 


0.178 


-0.194 0.7 


/T(«2) 


0.187 


-0.517 


0.014 


-0.117 


0.5 


0.040 


-0.762 


-0.201 


0.050 1.5 


/on«2) 


0.185 


-0.841 


-0.075 


-0.005 


0.4 


0.041 


-1.078 


-0.123 


0.068 1.2 


/t(«2) 


0.203 


-0.659 


-0.008 


-0.118 


0.3 


0.038 


-0.944 


-0.244 


0.073 1.5 


/f(«2) 


0.228 


-0.632 


0.017 


-0.143 


0.5 


0.049 


-0.931 


-0.245 


0.061 1.5 


/o''(«2) 


0.226 


-1.031 


-0.092 


-0.005 


0.4 


0.050 


-1.32 


-0.150 


0.083 1.2 


/#(a2) 


0.264 


-0.858 


-0.011 


-0.153 


0.3 


0.049 


-1.228 


-0.318 


0.095 1.5 


/:(a2) 


0.185 


-0.514 


0.014 


-0.116 


0.5 


0.039 


-0.757 


-0.199 


0.049 1.5 


/o''(«2) 


0.183 


-0.829 


-0.076 


-0.002 


0.4 


0.041 


-1.068 


-0.122 


0.069 1.2 


/^(«2) 


0.216 


-0.722 


-0.007 


-0.128 


0.3 


0.040 


-1.019 


-0.259 


0.076 1.4 


/TK) 


-0.141 


-0.775 


0.004 


0.161 


0.7 


-0.054 


-0.506 


0.621 


-0.326 5.2 


/on«4) 


-0.139 


-0.687 


0.170 


0.002 


1.5 


-0.061 


0.703 


0.323 


-0.209 2.9 


/?(a4) 


-0.167 


-0.895 


0.077 


0.143 


1.1 


-0.047 


-0.327 


0.698 


-0.394 4.9 


/^(«4) 


-0.173 


-0.947 


0.005 


0.196 


0.7 


-0.067 


-0.618 


0.759 


-0.398 5.2 


/o^(a4) 


-0.170 


-0.838 


0.209 


0.001 


1.5 


-0.075 


0.871 


0.392 


-0.254 2.9 


/#(a4) 


-0.217 


-1.165 


0.101 


0.187 


1.1 


-0.061 


-0.426 


0.909 


-0.513 4.9 


/:k) 


-0.140 


-0.770 


0.004 


0.159 


0.7 


-0.054 


-0.502 


0.616 


-0.323 5.2 


/o''(«4) 


-0.138 


-0.681 


0.170 


0.0005 


1.5 


-0.061 


0.710 


0.318 


-0.206 2.9 


/t(«4) 


-0.178 


-0.955 


0.083 


0.153 


1.1 


-0.050 


-0.349 


0.745 


-0.421 4.9 



Table C: Fit parameters for Eq. (A. 6) for the functions /"' defined in (A. 5). 5 is a measure 
of the quahty of the fit and defined in (A. 7). 



Split formfactors. As discussed in Sec. 4.2, the values of the Gegenbauer moments 
^1,2,4 are not very well known. In Sec. 4.4 and Tabs. A, B we have presented results only 
for our preferred choice of these parameters, i.e. 



af(l GeV) = 0.17 



■K,K,TI , 



1 GeV) = 0.115^ 



■n,K,ri I 



-IT,K,TI 



a, 

n,K,r] / 



1 GeV) = -0.015, 



af (2.2 GeV) = 0.135 , Oa "(2-2 GeV) = 0.080 , a4'"'''(2.2 GeV) = -0.0089; 

for set 4, the a^ are scaled up to /iir = 2.6 GeV. In order to allow the inclusion of 
future updates of these values, we split the formfactors into contributions from different 
Gegenbauer moments. We define^^ 



fiq') = r^(r) +ai(/iiR)r(r) +a2(m)r^(g') +a4(/iiR)r^(g^) 



fci4(^2\ 



(A.5) 



14 



Note that this sphtting is exact and vaUd for arbitrarily large a^ — there are no nonlinear terms in 
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ml (1-) 


mi (0+) 


ymax 


7r(r/) 
K 


5.322= 28.4 
5.412=29.3 


5.632=31.7 
5.722=32.7 


26.4 (22.8) 
23.8 



Table D: Masses of 1~ and O"*" resonances in the hu and hs channels. The 1~ masses 
are obtained from experiment and heavy-quark relations, the 0"'' masses from a potential 
model [40]. All numbers in units GeV2. 

where /"* contains twist-2 contributions from the asymptotic DA and also all higher-twist 



as,ai,a2)'^4 



m 



contributions not proportional to ai^2,4- The task is now to fit all functions / 
the interval < g2 < 14GeV2, to appropriate parametrizations. 

For /"*, which gives the dominant contribution to all formfactors, we use the same 
parametrisation as for the full formfactors. The results are collected in the lower halves 
of Tabs. A and B. Again, the fits are very good and best for the pion and set 2. 

The /"' turn out to be slowly varying functions of q^, which can be fitted by a poly- 
nomial of 3rd degree: 



p^q^) = a + h{q') + c{q'Y + d{q 



-,2\3 



(A.6) 



The measure of the quality of the fit has now to be defined in a slightly different way, as 
the /"' have zeros in the fit interval. We define the fit-quality 5 as 



100 



i:t\fit)-f\t)\ 

EJ/WI 



tG{0,i 



2' 



' 2 ' 



14} GeV2 



(A.7) 



i.e. as the average deviation of the fit from the true value, in per cent. The fit parameters 
are given in Tab. C As one can read off from the (5's, the fits are best for /"^, still good 
for f^'^ and worst for f""^. The limited quality of the fits for /""' is due to a change of sign 
of its derivative at the upper end of the fit interval, which cannot be reliably reproduced 
by a polynomial of 3rd degree. 

We would like to stress that none of the split-formfactor parametrizations must be 
used for q^ larger than 14GeV2. For calculating the full formfactors for arbitrary ai^2,4; 
the following procedure should be followed: 

• determine ai,2,4 at the scale ^\^ = m\ — ml; the scaling factors from /i = 1 GeV up 
to 2.2 GeV (i.e. m^ = 4.8 GeV) are (0.793,0.696,0.590) for (01,02,04); 

• choose set 2 (preferred) or set 4; 

• calculate /"* from the appropriate formula (A.l), (A. 2) or (A. 3), using the fit pa- 
rameters from Tab. A or B; 

• calculate f"-^'^'-^ from (A.6), using the fit parameters from Tab. C; 

• calculate the total formfactor from (A. 5); 

• extend the formfactor to the full kinematical regime by fitting it to (A.l), (A. 2) or 
(A.3). 
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B Distribution Amplitudes 

In this appendix we collect explicit expressions for all the DAs that enter the formfactors. 
These expressions are well known and have been taken from Ref. [16]. 

The key point is that, to leading order in QCD, DAs can be expressed as a partial wave 
expansion in terms of contributions of increasing conformal spin, the so-called conformal 
expansion. The coefficients of different partial waves renormalize multiplicatively to LO 
in QCD, but mix at NLO, the reason being that the symmetry underlying the confor- 
mal expansion, the conformal symmetry of massless QCD, is anomalous and broken by 
radiative corrections. 

The two-particle twist-2 amplitude (8) is expanded as 

0(n,/i) = 0,3H$^a„(/i)C^/2(C) (A.8) 

n>0 

with ( = 2u — 1 and Qq = 1 from normalization: 

1 
du(f){u, fi) = 1. 



The Cn (C) are Gegenbauer polynomials. The conformal spin of the term in Cj is 
j = n + 2. For the vr and rj one has a2n+i = due to G-parity, but af^ ~ {rris — rriq) for 
the K [33], which is one source of SU(3) breaking for the formfactors. 

As only the first few Gegenbauer moments a^ are known numerically, we truncate the 
series at n = 4; the values of the conformal spins included are listed in Tab. E, whereas the 
numerical values of the Oj are discussed in Sec. 4. The truncation is justified as long as the 
perturbative kernels T with which the DAs are convoluted are slowly varying functions 
of u, so that the rapidly oscillating Gegenbauers suppress terms with high n. In our case 
the T are nonsingular for all u, including the endpoints m = 0, 1, so the truncation of the 
series is justified. The term labeled 0as in (A.8) is the asymptotic DA which is reached 
for large scales // — *> oo; it is completely determined by perturbation theory and given by 

0as(M) = 6m(1-m); 

it is the same for all mesons. The Gegenbauer moments a„ become relevant at moderate 
scales and depend on the hadron in question. 

Let us now define the three-particle DAs. To twist-3 accuracy, there is only one: 

{0\u{x)af,ul5gGai3{vx)d{-x)\'JT~{p)) = 



rriu + md 



{PaPp^gvp - PaPvgt^p - PpPi^gua + ppPuga^d T(y, p ■ X) + . . . , (A.9) 



where the ellipses stand for Lorentz structures of twist-5 and higher and where we used 
the following short-hand notation for the integral defining the three particle DA: 

r(..p..)^/^,e-— .)r(„.„„.„.,. ,A.io) 
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Here a is the set of three momentum fractions ad {d quark), a„ [u quark) and Og (gluon). 
The integration measure is defined as 



/^^^/ 



daddaudag6{l — a^ — ad ~ cig)- 



There are also four three-particle DAs of twist-4, defined as 
(0|n(x)7^755fGc,/3(ra)d(-x)|7r^(p)) = 

= Pt,{PaXf3-ppXc,) — UmlA\\{v,p- x) + {P(39a,, - Pa9pf,)UmlA±{v,p ■ xOA-H) 

{0\u{x)-f^igGapivx)d{-x)\7r~{p)) = 

= PtiiPaXp - Ppx^) — UmlVii{v,p- x) + {ppgi^ - Pa9if,)fnmlV±{v,p- x)^A.12) 



g^^, is defined as 



V = 9t^u [Pt^x^ + p^Xf,). 

^ p ■ X 



To next-to- leading conformal spin (j = 7/2,9/2), the twist-3 three-particle distribution 
amplitude T is given by^^ 

T(a„, ad, ag) = 360?73a„adag{l + uj3-{7ag - 3)} . 

The two-particle twist-3 distribution amplitudes (j)p and 0o- in Eqs. (9) and (10) depend 
on T through the equations of motions [16],^^ which implies that their coefficients are not 
independent from each other. The expansion up to NNL order {j = 3/2, 7/2, 9/2) reads^^ 

Mu) = 1 + {30r/3 - Iplycl^'iO + {-3r/3^3 - ^P^ - ^P^a2}Cr(C), 

1 7 3- 

(j)^{u) = 6m(1 -u){1 + {5?73 - -r/3^3 - —pi - -p^agjcf ^(C) . 

The two-particle twist-4 corrections g.^ and A in Eq. (8) are given to NNL conformal spin 

^^In the literature the notation /s^ = fTrVs is also widely used. 

^^An explicit expression for (j)p in terms of T is given in Ref. [41], Eq. (16). 

^^At first glance it seems that (fip is taken to a higher order in conformal expansion than if)a-, but as 
discussed in the first reference of [16], (f>p and (f>a- arc not pure spin projections, which means that the 
coefficients of a given Gegenbauer polynomial contain contributions from different partial waves. 
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(j = 1,3,5) by 



18 



18 „„ 20 .^1/2,^, , , 9 „ .^1/2, 



^.(n) = 1 + {1 + -a2 + 60r/3 + yr/4}C2''lC) + {-^^2 - GvsiOsjC.^^C), 
w X _ f 16 24 20 

A(m) = 6mm < — + — 02 + 207/3 + — ^74 



18 \ 

— 02 + 21r/4CJ4 1 {2u^(10 - 15n + 6^^) Inn + 2u-'^(10 - I5u + 6u^) \nu 

+uu{2 + 13uu)} . 
Finally the three-particle twist-4 DAs are to NL spin (j = 3, 5) given by 
^ll(a) = 120a„adag(aio(ad- «„)}, 
V||(a) = 120auadag{voo + Vio{3ag - 1)}, 
A±{a) = 30a^(a„ - ad)[hoo + /loictg + hio{5ag - 3)/2}, 

V±{a) = - 30aJ{/ioo«9 + /ioiK^g - 6a„aJ + /iioK^g " 3/2(a^ + a^)]} , 

where a = 1 — a and the aij, Vij and /ijj are related to hadronic matrix elements 774, u^ 
and a2 as 



aio = 


= fr/4^4- 


- 20«2, 


WlO = 


= Y^4UJ4 


Woo = 


= -|r/4 


/ioi = 


= |r/4u;4- 


■ 20«2 , 


/ilO = 


= IVi^A + ^02 , 


^^00 = 


= -|r/4 



Taking everything together, we have 7 hadronic parameters {q} = {ai, a2, 04, 773, ci;3, 774, ci;4} 
which parametrize all DAs to twist-4 and NLO in conformal spin. The q are scale- 
dependent and are usually given at the scale 1 GeV. To LO in QCD, they do not mix 
under renormalisation, so that the scaling up to /xir = ^^m^^^m^ is given by 

Q(/iiR) = L>«/*Q(lGeV), 

with L = as(/iiR)/as(l GeV), j3o = 11 — 2/3Nf. The one-loop anomalous dimensions 
7cj are given in Tab. F. Note that the anomalous dimension increases with increasing 
conformal spin, 7 ~ log j, which implies that the truncation of the conformal expansion 
becomes the better the high the scale. The numerical values for all these parameters at 
the scale /i = 1 GeV are collected in Tab. G, taken from Ref. [16]. 

""^^Note that, contrary to appearances, the contributions oi Qt^ and A to (8) do not vanish for zero meson 
mass: 7% imphcitly contains a factor l/rnj and survives in the Hmit m^ -^ 0. 
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tree 


0(a.) 


twist 


2 


3 


4 


2 


3 


4 


x-particle 


2 


2 


3 


2 


3 


2 


2 


3 


2 


3 


JL 
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Jnnl 


2 
4 
6 


3 

I 

2 


V 
2 


1 
3 

5 


3 

5 


2 

4 
6 


3 

? 
2 


- 


- 


- 



Table E: Overview of the contributions included in the calculations. For the K we also 
include conformal spin j = 3 for twist-2 which explicitly parametrizes SU(3) flavor break- 
ing. 



7a,i 


Irr, 


'7a;3 


7r)4 


li^i 


Cf(i- 


2 

{n + l){n + 2) 


n+l \ 


fc.c. 




3^^ 





Table F: One- loop anomalous dimensions of hadronic parameters in DAs. 





IT 


K 


V 


m 


0.015 


0.015 


0.013 


tUs 


-3 


-3 


-3 


^4 


10 


0.6 


0.5 


UJ^ 


0.2 


0.2 


0.2 



Table G: Input parameters for twist-3 and 4 DAs, calculated from QCD sum rules. The 
accuracy is about 50%. Renormalization scale is 1 GeV. 



C Spectral Densities for / 



+ 



The total spectral density of 11+ is obtained as sum of all the contributions listed below, 
i.e. 

Pn+ = PT2 + PT3 + Pa + Pp + Pt4 + Pt4- 

Pt2 is the contribution from the twist-2 DA, pt3 from the twist-3 three-particle DA, p^(^p) 
from the twist-3 two-particle DA (p^ip) and Pt4 from the two(three)-particle DAs of 
twist-4. There is also one constant term, T4c, which is due to twist-4 corrections that 
cannot be expressed via a dispersion relation, so that the total Borel-transformed 11+ is 
given as 



J ml 
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PT2 



(g2 - s) 



^4 „2 , 1 n „2 „3 , 4 



j{mt - q^) {mi - s) (15 04 (42 mg + g« + 10 g% + 20 q^ s^ + 10 g^ s^ + 



- 84 mg (q2 + s) + 28 m^ (2q2 + s) [q^ + 2 s) - 14mg (g^ + s) (g"* + 4g2 ^ _^ ^2^^ ^ ^^2 _ ^-)2 
x{&a2{bml + q^ + iq^ s + s^ -bml {q^ + s)) + {q^ - s) (oq (g^ - s) + 3 ai (-2 mg + q^ + s)))) 



+ a4 ^^°/-";3^ (mg-.)(-2mgg^ + 2g4 + mg(4 + 7r^),s-(l + ^^)g^s-3,s^) + (mg-a 



2 _ „^3 



s(g2 — s) 



X (s (s - 771^) log(l - -^) + s log( — ) (-2 s + (m^ - s) log( — )) + 2 (777^ - s) (2 777^ - 5 s 

777i 777f 777f 



+ 2 s log(-3)) log(^ -l) + 2s{s-ml) log(^ - 1)" - 2 (777^ - s) log(l - ^) {ml 



777.^ 771^ 



+ S + S log( — ) - 2 s log(— - 1))) + 2 (7772 - q^) (ml - s) (Li2(- 



■ 9 9 ' 

q'^ — m^ 



Li,(l-^) 



+ 4Li2(l - 4))) + ^ll^iiml - s) (6g2 (6,77^ - 6777^^2 ^ q'') + [-2m\ (95 + 9^^) 

+ 27 777^ (6 + TT^) q^ + {3-9 n^) q'^) s + (777^ (158 + 9 tt^) - 3 (50 + 3 tt^) q^)s^ - 13 s^) + 3 (8 (777^ - s) 



X {q' - mi) (2 m^ - q' - s) s log(^) + 3 {mi - q') {ml - s) (2 777,2 _ ^2 _ ^) ^ j^g^^ _ ji^^ 

mt mf 



s 
'ml 



+ s log(^) (12 777^ {q^ + s) + 14g2 s (g2 + s) - 2777^ (4g4 + 19g2 s + 7s2) _- 3 {ml - q^) {ml - s) 



X (2 7772 ~ g2 _ ^) log(^)) - 2 (777,2 _ ^2) (^ __ „^2) (_2 7772 + g2 + s) (6 ^2 _ 23 5 + 6 s log(^)) 
mi mi 

X log(^ - 1) + 6 (7772 - g2) s (s _ ^2) (_2 ,772 + 92 + s) log(^ - 1)' + 6 (7772 - g2) 



X (7772 - s) {2ml -q^-s) log(l - ^) (m2 -s + s log(^) - 2 s log(^ - 1))) - 18 (7772 - q^) 

m^ 777j, m^ 

x{ml~s){2ml-q'-s){U,{^^)+Lh{l-'^) + AU,{l-^))) 

g^ — 777^ s m^ 



a2 /tt "^b ^f_2 „\i o^ ^_2 „2^„2^Qn™4 ic;™2„2 , „4^ , ^c; ^6 m qq , 70 „2^ 
4s(q2_s)5 



-((777^ - s) (-24 (777^ -q^)q^ (30 777^ - 15 777^ q^ + q^) + (5 777^ (1183 + 72 7r2) 



20 777^ (407 + 36 7r2) q^ + 12 {b - Q i:'^) q^ + 216 ml {II + 2 n^) q"^) s-{5mf, (1525 + 72 tt^) 

16 777^ (575 + 36 7r2) q^ + 36 (73 + 6 n^) q^ ) s^ + {ml (2083 + 72 n^) - 8 (260 + 9 tt^) q^) s^ - 61 s^) 
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+ 12 (25 {ml ~ g2) (mg ~ s) s (5 m^ + g^ + 3 q^ s + .s^ - 5 mg (g^ + .5)) log(-^) - 6 (mg - q"^) 

ml 

1 2 
X [ml - s) s (5m^ + g4 + 3 g2 s + .,2 _ 5,^2 ^^2 ^ ..^^ j^g^^ _ 9 ^ _^ _^ log(^) (-60 ml [q^ + s) 

m^ m^ 

+ 37q^s [q^ + iq^ s + s^) + iQml{?,q^ + ^q^ s + is^) - ml (25 q^ + 237 q'^s + 2Qlq^ s^ + 37 s^) 
+ 6 {ml- q^) {ml - s) {bml + q^ + iq^ s + s^ - bml {q^ + s)) log{-^j) + 2 (mg - q^) {ml - s) 

x{5mt + q'^ + 3q^s + s^ -5ml{q^ + s)){l2ml-55s+12s log(^)) log(^ - 1) + 12 s 

X (mg - q^) {s - ml) {bmt + q"^ + Sq^ s + s^ - 5 ml {q^ + s)) log(^ - 1) - 12 (m^ - q^) 

ml 

2 

X {ml - s) (5 m^ + g^ + 3 g^ s + s^ - 5 m^ (g^ + s)) log(l - -^) {ml - s + s log(^) -2s 

m^ m^ 

X log(^ - 1))) + 144 {ml - q^) {ml - s) (5 m^ + g"* + 3 g^ .s + s^ - 5 ml {q^ + s)) 
m.f 



X (Li2(^^^) + Li2(l - -^) + 4Li2(l - — ))) 
q^^ — m^ s ml 

+ ^±ljl^^{-{{ml - s) (30 {ml - q^) q^ (1260m« - 1890mgg2 _ lOSmg^e + 2g« 
103(52 — s) 

+ 840 mlq^)- (21 m^° (23207 + 900 tt^) - 63 mg (18827 + 900 tt^) 9^ _ 700 ml (439 + 45 tt^) g^ 
+ 150 ml (176 + 45 tt^) q^ + 30 (19 - 15 tt^) q^° + 175 mg (5603 + 360 tt^) g^) 5 + 54 {ml 
X (13051 + 450 TT^) - 112 mg (22157 + 900 tt^) g^ „ goO ml (616 + 45 tt^) g*^ + 150 (283 + 30 tt^) g« 
+ 525m^ (3523 + 180 tt^) q^) s^{-lml (119363 + 3600 tt^) + 756 mf (2131 + 75 tt^) q^ + 200 (848 
+ 45 TT^) q^ - 10125 ml (89 + ATT^)q^) s^ + {7 ml (34967 + 900 tt^) - 6 ml (57989 + 1800 tt^) g^ 
+ 125 (1075 + 36 TT^) g4) s4 ^ („2 mg (10553 + 225 tt^) + (21461 + 450 tt^) g^) ^5 _^ -^g-^ ^6^^ ^ gg 
X (91 {ml - (7^) {ml - s) s (42 m^ + g^ + 10 gS + 20 q^ s^ + 10 g^ s^ + s^ - 84 m^ (g^ + s) 

„2 

+ 28 m^ (2 g2 + s) (g2 + 2 s) - 14 mg (g^ + s) (g^ + 4 g^ .s + s^)) log(-^) - 15 {ml - q^) {ml - s) s 

ml 

X (42 TO^ + gS + lOgS + 20g^s2 + 10 g^ s^ + s^ ~ 84 m^) (g^ + s) + 28 m^ (2g2 + s) (g^ + 2 s) - 14mg 

2 2 

X (g2 + s) (g4 + 4 g2 s + s^)) log(l - 1^) + s log(4^) (-1260 mj" (g^ + s) + 630 ?7ig (5 q^ 

mt mt 
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-12 



q^s + bs^)~ 210 ml {q^ + s) (13 q"^ - 



' ml (91 q^° 



3 g'' + 48 q^ s + 13 s^) + 121 q^ s {q^ + 10 q^ s + 20 g'* s^ 
+ 10 q'^ s^ + s^) + 105 ml (9 g^ + 82 g*^ s + 160 g^ s^ + 82 g^ ^3 ^ g ^4^ _ , 
+ 2305 q^s + 9400 q^ s^ + 9700 g" s^ + 2575 g^ ^4 _^ ^21 s^) + 15 (mg - g^) (^2 _ ^-j 
X (42 m^ + g* + 10 g*^ s + 20g''s2 _,_ ^q^^ ^3 _^ ^4 _ g^^^^e (-^2 ^ ,,-) _^ 28 m^ (2?^ + 5) (^2 ^ 2 s) 
- 14 mg (g2 + s) {q^ + Aq\s + s^)) log(^)) + 4 (m^ - g^) („^2 _ ^) (43 m^ 
+ 20g'*.s2 + 10 9^ s^ + s'^ ~ 84 mg (g^ + 5-) _^ 28 m^ {2q^ + s) (g^ + 2 s) - 14mg {q^ + s) 

2)) (15mg - 83 s + 15 s log(^)) log(^ - 1) + S 

mjj 771^ 



^tg + 9^ + 10gS 



; (9^ + 492 s- 



X (42 ml + q^ + 10 g^ s + 20 <?"* s^ + 10 q'^ s^ + s"* 



imt{q^ + s) 



30 {ml -q^) sis -ml) 
28m^(2(72 + s)(g2 + 2s) 



X 



^2 (g2 + s) (g4 + 4 92 5 + s^)) log(^ - 1)' - 30 {ml - q^) {ml - s) (42 mf + g^ 

m^ 

m 

2 
(9' + s) (g' + 4 g" s + s')) log(l - ^) (77^2 - s + s log(^) - 2 s log(^ - 1 

777^ 777^ 777^ 



-14 777^(^2 

+ 10g^s + 20g''s2 + 10g2s3- 



-s"- 



84mt{q^ 



{q^ + s) + 28 777^ (2 (j2 + s) (g2 + 2 s) - 



14 777^ 



1))) 



+ 900s (777^ - q^) {ml - s) (42 777^ + g*^ + 10 g*^ s + 20 g'' s^ + 10 q^ s^ + s* - 84 ml {q^ + s) 

2 2 

+ 28 7^^ (2g2 + s) {q^ + 2 s) - 14 7ng (g^ + s){q^ + Aq^s + s^)) (Li2(^^ 2) + ^12(1 - — ) 

g^ — 777^ s 

+ 4Li2(l--^))} 

— 15 7^3 ^T 2 2\2 / 2 \/r,/2 \/r 2,o2,o 

PT3=— T ^("ife-g) (777b -s) (2 (g -s)(-5777fc + 2g +3, 

(r - s) 



\ I f'7 i /3 22,4 o2 
)) + (7 777(j — 6 771^ g +q — 8 777(,S 



f4(72s + 2s2)cJ3) 
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Pp 



ao^i^^ {3mls -3q^ s) brjs^: 



6 (g^ — s) s 



^(18 ml s - 36 mf q"^ s-3q^s + 21 ml q^ s - 18 mf s^ 

2 — s) s 



30mU^ s^ - 12q\s^ + 3mls'~3q^ s') + a,{ ''° ^" , ((g^ + s) {s -3ml) 

6 (q^ — s) s 



q_ 

m 



2 2 

2 
b 



-3{ml-4q^)s\og{\} -18s{3ml~3q^ + s)log{l-^) (ml - i q^) s \og{^) 






,s-q 



X (log(l - ^) - log( 2^)) + 6 log(l - ^) (m^ <7^ + 3 m^ s - g^s + s ((2 ml - 3 q') 



log(-^) + g2 log( — )(m2 + g2 + ,) log(^^) - (mg ^ 2 q^) log(— - 1))) + 3 (m^ 



S 2 

mt 



Xlog(^) -log(^-^) s(5g2 + s) +log(^)(2m2g2+^2^_352^_4^2^jQg^^_-^)) 



s 

TO? 



S 

ml 



2 2 

2 log(^^^) (2 g2 s _ ,7^2 ^^2 + s) _ s (2 mg - 3 g2 + s) log( Aj) + s (3 mg - 2 q^ + 2 s) 



X log(— - 1)) + 2 log(— -l){2q's- ml (q'+As) + s (s log(^) - (to^ + 2 s) log(— - 1)))) 
TTif m^ TO,f mf 



6 s {ml Li2(-^-^ 2 ) + (2 ?7i2 - 3 9^ + s) Li2(l - -4) + to^ Li2(l ^) - (mg + g^ + s) 

(7^ — ?7l^ toJ S 



X Li2(l - :^) - (5?^ + s)Li2(^^^i— 2 ) - K + 9' + •^) Li2(^— ) - ("ifc + g' + s) 



mg -r 



g' -s 



xLi2(^47^) + (2"^b-3g^ + .s)Li2(^^^^) + 2(TO2-g2)Li2(i-4) + (™2_4^2) 



X Li2(l - ^)))6 (g2 _ 5)' .s + -i!^l^(-3 TO^ {q' - 12 to^) g* - 6 to« (^^ - 47) s - 18 to* (24 



s 
? 



(9^ "5)5 



+ TT^) g2 s 4. (3 4. 2 TT^) g«^ s + 3 7«g (8 TT^ - 17) g" s + 6 TO^ (tt^ - 60) 3"^ + 3 ml (205 + 6 tt^) g^ s 



(3 + 26 TT^) g4 s2 _^ gg „^2 ^3 _ (--^47 ^ 3 TT^) q^ s^ + (-9 + 2 TT^) s* + 3 (to^ - 4 g^) (g2 ^ s)^ s log(-^) 



nit 



2 2 2 2 

+ 6 {ml - 4q2) (q2 - sf s log(^) (log(l - ^) - log(^^)) + 6 s log(-^) (15 ml - 15 to* (2 q^ + s) 

TO,^ 771^ 771^ ?Tt^ 

-3g^s(g2+4s) + m^(9g'' + 30g2s + 6s2) + (6mj^-6TO^s-3g2(g''-4g2s + s2)+TOg(-4g*-4g2s 
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mi 



+ 2 s^)) log(l - ^) - {q^ - sf (2 m2 - 3 g2 + .,) log(^^) + (-6 mg + 6 m^ g^ + 6 m^ s + s (-5 q^ 



2q^s + s^) ) log(^ - 1)) + 3 (-(s (6 m^ - {q^ - s)' {3q^-s)-6 m^ {q^ + s) + 3 ml {q^ + s^)) 



2 2 

log(l - ^) ) + m^ s (6 m* + g4 + 4 g^s + 5^6 mi (q' + s)) log(^) - [q' - sf log(^^) 
mf mi mi 



X (2 (-2 g2 s + 771^ (g2 + s)) + s (5 q^ + s) log( ^))2 (6 mg (2 g^ + 9 s) _ 6 m^ (2 g* + 18 g^ s + 9 52^ 

- 2 g2 s (g4 ^ 7 g2 5 _^ ^3 ^2^ _^ ^^2 (-^6 + 44 g4 ^ + 89 g2 52 ^ ^g s^)(-3 mg + 2 g^ _ 2 s) (g^ - sf s 

2 

X log(^^^)) log(4j - 1) - 2 s (-6 m? - 6 mf (g^ - s) + 2 s (-5 g'' - 2 g^s + 5^)771? (13 g^ + 4 g^ s 
™6 7772 

+ s^)) log(^ - 1)^ + log(^) (-12 777^ s^ _ 3 q2 .. (^2 ^ 3 ..^2 _^ ^2 (2 ^6 ^ 21 g4 ^, + ^2 g^ s^ + 13 s^) 
r?7,f m,t 



- 4g2 s (6 777^ + g4 + 4g2 s + s^ - 6 777^ (g^ + s)) log(^ - 1)) + 2 log(l - -^) (6777^ (g^ + 2 s) 

- g2 s (g* - 8 g2 s + s^) _ 6 777^ (g-* + 2 g2 s + 2 s2) + 777^ (g^ + g^s + g^ s^ + 3 ^3) ^ ^ ^^2 ^g ^4 _^ ^4 



s 
777; 



+ 4 g2 s + s2 - 6 777,^ (g2 + s)) log(^) + (g^ - sf {ml + q^ + s) log(^^) + (-12 777^ g^ + 777^ (11 g^ 



+ 14 g2 s - s^) + 2 g2 (g" - 8 g2 s + s')) log(^ - 1)))) + 6 s (777^ (6 7n^ + g" + 4 g2 s + s' - 6 mf, 



X (g2 + s)) Li2( ,^ J + (g' - sf (2 777^ - 3 g^ + s) Li2(l - ^) + 777^ (6 777^ + g^ + 4 g^ s + ^2 
g^ — 777^ 7)7^ 



-67^? (^2 _^ ^)) Li2(l - ^)(g2 - s)' (777? + g2 + s) Li2(l - ^) - (g^ - ,s)' (5 g^ + s) Li2(^^^4) 

s mf^ — g^ 



2 , 9 9 . s^. .TO^ -s^ , -7 x2 , 9 9 . s^. (g^ -'•"ps 
q^ — s ' " ... - . j.j^2 ^qi _ g^ 



(g2 - sf {ml + q^ + s) l^U^^^^-^) - {q' - sf {ml + q' + s) Li2( 



) + {q'-sf{2ml 



- 3g2 + s) U2C-^ 1-) + 2 (777^ - g2) (6777^ + g^ + 4g2 s + s^ - 6777^ (g^ + s)) Li2(l - A) + ("^? 

s — q'^ mi 



4g2)(g2-srLi2(l--))} 
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P- = ., . ^-. ("" ('^' - ^) (-("^' (3 g2 + s)) + g2 (q2 ^ 3 ..)) _ 30 ^3 (10 ^8 + 10 ml {q^ - s) 

- 30 ml g2 (^2 ^ .,) ^ ^2 (13 ^6 ^ 45 ^4 ,, _^ 31 ^2 s2 ^ ^^3) _ g2 (^6 ^ 11 ^4 ^ _^ 15 ^2 ^2 _^ 3 ^3))) 
2 

6 (q2 _ s) s 



+ 4 (-15 + TT^) 9^) s^ + (3 + 2 TT^) 3-' + 6s (-4 g^ s + mf (3 g^ + s)) log(l - ^) + 6 s log(-^) 

rn^ m^ 

X (3 (-2 q^s + ml {q^ + s)) + {q^ - s) {-{{ml + q') log(l - ^)) - {ml + q' + s) log(^^i^) 

mf q'^ — s 



2 2 

,2 ,1 „2^ i„„/™b ~ 9 ^ , ^_2 , „2 



^ -' -'2 „\ „ /_2 , „2 -■ ' ■* 



- {mi - 4g^) log(^-^) + {mi + g^) log(— - 1))) + 3 (-2 (g^ - s) s {mi + q' + s) log( — ) 



s — q^ 



mt 



+ 2{q'- s) log(^^^^) {-{^l {q^ + 2s)) + s {ml + q' + s) logC^-f)) + {q' - s) s {{ml 



, mg — s , 
g2 — s 



9 s , , f7if — g" , ,9 9 . , , ( Tif — g ) s , „ , 9 9 , , . \m,t ~ q ) s . 
+ g^ + g log " ^ - mg + g^ + ,s log ^ / ^ - 2 mg + g^ + ,s log ^ ^^ 2^ 

s — g^ 771^ {s — g^j m^ (s — g'^j 



,s-g- 



X log(:-^) + 2 (2m2 - 3g2 + s)\og{'—^) + 2 log(^^^^) (-3 77^^ + Gg^ + (2 777^ 



mt 



s-q- 



. 3g2 + ,s) log(-^))) - 2 (4777^ (g2 + ,s) + 2g2.s (g^ + 3.s) - 27772g2 (g2 + 7s) + {q^-s). 



X (5 7^2 _ 5 g2 + 3 s) log(^^)) log(^ - 1) + 2 s (-2 s (g^ + s) + Tn^ (3 g^ + s)) log(^ - 1)' 

TTlf 777i ?77f 



+ log( — ) {q' s{3q' + s)- mi {q' - s) (2 g^ + 5 5) - 2 (g^ - s) s {mi + g^ + s) (log( " I ) 



,ml-q^ 



s-q- 



-\og{^<-f}^)-log{'-^)) + 2s{-5mlq^ + q^~3mls + 6q's + s')log{^-m 
m^i {s — q'^) mi mi 

o2 
+ 6 log(l - -^) (-(g2 (g2 - 7 s) s) + 2 777^ (g2 + 2 s) - ml (g^ + 9 g^ s + 2 s^) + s ((4 g^ s - ml (3 g^ 
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mt 



+ s)) log(^) + (2 7^g (g2 + s) g2 (g^ + 3 s)) log(-^) + (g^ - s) {-{{ml + g^ + .) (log(::fi^-f ) 



,7772 -g2 



s-g^ 



S 
ml 



K\ 'b ))-K-4g2)log(^^))-2(g2(3g2-7,s) + 777|(g2 

777^ (S - g^j 777^ 

+ 6 s ((4g2 s - 7772 (3 g2 + ,s)) Li2(l - ^) - K - 4g2) (g^ _ s) (2Li2(^^|-^) + Li2(::^-:) 



777J -g- 



g--s 
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-Li,( \^"^^ 'l ))-(3ml{3q' + s)-{q' + s)i5q' + s))U2{l-^)))+ J'!'''\, (18mgg^ 



'ml [q^ -s)' 



T-b 6 (g2 - s) s 



X (40 ml - 40 ml q^ + 9'') + (5 ml (-359 + 24 tt^) - 30 mf (127 + 16 tt^) g^ + 30 ml (129 + 20 tt^) 9" 
+ 24 mg (40 - 11 TT^) gS _^ 6 (-3 + 2 TT^) g8) s + 2 (5 ml (295 - 36 tt^) + 30 m^ (129 + 16 tt^) g^ _ 6 mg 
X (609 + 68 TT^) g4 + 6 (-27 + 20 tt^) q^) s^ + 6 (5 m^ (-25 + 8 tt^) - 19 ml (39 + 4 tt^) g^ + 9 (55 
+ 4 TT^) g4) s3 _ 4 (-^^2 (71 + 6 TT^) - 66 g2) s^ + (-37 + 12 tt^) s'^ + 36 s (10 mg - 10 mg s - 10 ?7i^ g^ 



X (g^ + 2 s) - 4 g^ s (9* + 3 g^ s + s^) + m^ (3 g'' + 25 q* s + 21 q^ s^ + s'^)) log(l - -^) + 12 s 



TOt 



X log(-^) (-(m? (5 m? - 2g2) (5^4 _^ 20 m? g^ - 28 g'')) + (40 m? + 135 mtg^ _ 249 m? g'' + 47 g^) s 
ml 

- 3 (5 m^ + 18 ml q^ - 32 q"^) s^ + 9 {ml - 2 g^) s^ - 3 (10 mf - 20 mf (g^ + s) + 10 ml (g^ + s) 



X (2 g2 + s) + g2 (^6 + 7 g4 , + 3 ^2 ^2 _ ^3) _ ^2 (9 g6 ^ 23 g^ s + 7 g^ s^ + s^)) log(l - -^) 



mf 



+ 3 (g2 - ,s)' (-((m^ + g2 + s) log(^^^)) - {ml -4q') log(^^^^4-)) + 3 (10 m« - 20 m^ (g^ + s) 

"^ " s — q'' 



q^ — s 



+ 10 ml {q' + s) (2 g^ + 5) + g^ (g" + 7 g* s + 3 g^ s^ - s'^) - m^ (9 g" + 23 g^ s + 7 g^ s^ + s'^)) 



log(^ - 1)) + 36 log(l - -^) (20 ml g^ - 10 ml g^ (3 g^ + 5 s) + 2 7«^ (g^ + 2 s) (6 g* + 13 g 



mr 



mr 



■ s^) -q^s{q^ + q^s- 15 g^ s^ - 7 s^) - ?7ig (g^ + 7 gS + 35 q^ s^ + 35 q^s^ + 2 s^) + s ((-10 mg 



+ 10 mg s + 10 mt q' (g^ + 2 s) + 4 g^ s (g* + 3 g^ s + s^) - mi (3 g'^ + 25 g* s + 21 g^ s^ + s^)) log(-i^) 



+ (10 ml (2 g^ + s) - 10 m^ (3 g'^ + 5 g^ s + s^) + 2 m^ (6 g" + 24 g* s + 14 g^ s^ + s"^) - g^ (g*' + 11 g* s 



. ; - W - ^)' {-{{ml +q' + s) (log(I^^) - log(i!^;L^ 
^^ s — q'^ m-l {s — g^j 



+ 15 q' s' + 3 s')) log( — ) + (g^ - s)-- (-((m^ + g^ + s) (log(^^^) - log( 



))) - {ml 



- 4 g2) log(^-^)) - 2 (10 m? (2 q^ + s) - 10 mf (3 g* + 5 g^ s + s^) + g^ (3 g^ - 23 g-* s - 3 g^ s^ _ 7 s^) 
mt 



+ ml (11 g*^ + 51 gS + 25 g^ s^ + 3 ^3)) log(^ - 1))) + 6 (-6 (g^ - sf s {ml + q^ + s) log(^) 



S_,2 
2 - 



1 

q^ — s 



+ 6 (g2 - s)' log(^^^) (-(m^ (g^ + 2 ,s)) + s {ml + g^ + ,s) log(^^)) + 3 (g^ - sf s {{ml + g^ + s) 
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xiog(^^^)' - {ml + q' + s) iog[^^4r^S -n< + e + s) iogfi4r-^0 i°g(^) 

s — q"^ m^{s — q^) m^[s — q"^) m^ 

2^22 2 

+ 2i2ml-3q^ + s)logC-^) +2 1og(^^^^^— ^)(-3m2 + 6g2 + (2m2-3g2 + s)log(^^))) 

m^ s — q'^ m^ 

+ 2 {ml (-120 g^ + 170 s) + 20 ml (9 g"* + 24 q^ s - 10 s^) - 6 ml (12 g^ + 148 gS + 123 g^ s^ - 3 s^) 

- 2 g2 s (3 gS + 98 g^ s + 195 g^ s2 ^ 9 .,3) + 2 mg g^ (3 ^6 ^ ^^40 ^4 ., ^ 53^ ^2 .,2 _^ 216 s^) + 3 ^ (., _ q^f 

2 

X (5 m2 - 5 g2 + 3 s) log(f— ^)) log(4j - 1) - 6 s (10 mg - 10 m^ (4 g^ + 3 s) + 2 s (g^ + s) 

ml ml 

X (11 g4 - 2 g2 5 + s^) + 10 m^ (5 g"* + 8 g2 5 ^ 2 s^) - m? (23 g^ + 65 9* s + 41 g^ 52 _^ ^sxx iog(4j - 1) 

ml 

+ log(^) (-6 m? 9^ + g^ (120 ml - 157 m? g^ + 9 ^4^ ^ _|_ (^gQ ^6 _^ gg ^4 ^2 _ 343 m? g^ + 205 g^) s^ 
ml 

2 2 

- 3 (10 ml + 43 ml g2 - 121 g^) 5^ + 3 (5 mg + g^) 5^ _ 6 (g^ - s)' 5 {ml + g^ + 5) (log( "''' ' ^ ) 

s — g^ 

- log( ^"!'~^'^' ) - log(^^)) + 6 s (g« + 24 gS + 30 gS^ + 4 g^ s^ + s^ - 20 mg (2 g^ + s) 

771^ (s - g2) m;^ 

+ 20 mf (3 g^ + 5 g2 s + s^) - m? (25 g^ + 93 g^ s + 59 g^ s^ _^ 3 s^)) log(4j - 1))) + 36 s ((-10 m? 

m.l 

+ 10 ml s + 10 TOfc g^ (g^ + 2 s) + 4 g^ s (g" + 3 g^ s + s^) - m^ (3 gS _^ 25 g" s + 21 g^ s^ + s^)) Li,(l - -^) 

ml 

-(mg-4g^)(g^-.f(2Li,(^4^)+Li,(4^)-Li2( ''irr1 ))-(10"^g-5g^-16g^^ 

m^ — g q'^ — s 771^ (g^ — sj 

- 54g'^s2 - 4g2 s^ - s'' + lOTTig (4g2 + s) - IOttiI (7g'' + 12g2 s + 2s2) + mg (29g^ + 115g^s + 83g2 s^ 

+ 3s3))Li2(l-4))} 

r 2 

^2p ^ UrUbH^ ^ ^27 ^2 (2880 mj" - 25 m^ (215 g^ + 229 s) + 20 ml (149 g"* + 447 q^ s + 184 s^) 
120 (g2 — s) 

- 10 ml (38 g^ + 352 g"* s + 427 g^ s^ _^ 83 s^) + 5 g^ (g« + 17 gS + 42 q^ s^ + 22 g^ s^ + 2 s^) 

+ ml (-76 g^ + 84 gS + 584 g^ s^ + 584 g^ s^ + 24 s'^)) + 5 (16128 r/s ml° W3 - 2520 773 mf (15 g^ 
+ 13s)u;3 + 8g2 ((3 + 90r/3 + IO774) (g^ - sf {q^ + 5q^s + 2 s^) -9^ (g^ + 17gS + 42 g^ s^ 
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+22 q^s^ + 2 s^) tjg) + 3 mg (20 q^ s^ (3 - 16 774 + 144 (773 + 2 773 W3)) + q^ s (69 + 320 774 + 240 773 (6 + 19 LU3) 
~ 420 774 UJ4) + 3q^s^{-9 + 240 7/3 (-2 + 3 W3) - 140 7^4 W4) + 5 g* (-15 - 16 7/4 + 144 773 (-3 + W3) + 84 774 W4) 
+ s^ (-27 + 80 774 + 240 773 (-3 + ^^3) + 420 774 W4)) + 5 77ig (s^ (-33 + 352 774 + 144 7^3 (-8 + 33 W3) + 756 774 W4) 
+ q'^ (-33 + 352 774 + 144 773 (-8 + 47 W3) + 756 774 W4) + 2 q^s {33 + 576 7^3 (2 + 11 ^3) - 4 774 (88 + 189 W4))) 
+ 3 7n^ (-15 g2 s2 (1 ^ 48 (,^g _^ 16 773 u;3) - 12 774 (4 + 7u;4)) - 3 g"* s (53 + 1680 (773 + 3 773 W3) - 20 774 (4 + 21 1^4)) 
+ g*^ (111 + 240?73 (15 -19^3) -20774(20 + 63^4)) + 5^(240 773 (9 -IOW3)- 7 (-9 + 20 774 (4 + 9tJ4))))) 

9 9 ? 

Ill ~~ 5 JTl — Q 

+ 1080 77^2 {ml - q^) {ml - s) (2 77^^ - g^ _ ^^ (^2 _ ^) (g ^^ _ 35 ^^ ^^^ (_ j^g^^^ ^ ^ log( '' , ))) 

g^ — s s — q^ 



r 2 

^3p _ -Umb^J.^ ^_2 _„2w8n^. ^„2 _ ,^,'_a^4_L„4 _ q „2 „ _ « ^2 ^ _2 ^„2 _l 1 c; ,^^ _l 07„. r/ic;^6 ^ „6 _ 4„4 

24 (g2 _ s)^ 



p^P = ^;7^ 6 ("^b - 9') (80^4 (g^ - s) (-8 771^ + g4 - 3g2s _ 6s' + m^ (g^ + I5s)) + 27 as (45 77if + g^ - 4gS 



- 30 g2 s' - 12 s^ - 5 77^^ (7 g2 + 20 s) + 77ig (g^ + 68 g^s + 66 s'))) 

T4e = djLj:!il!^e^ (12 a2 + 5 (-25 + 48 7,3 (-10 + u;^))) 
160(7772 — q"^) 
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